AD-A018  281 


EVALUATION  OF  H I GH- RESOLUT I ON  EARTH  RESISTIVITY 
MEASUREMENT  TECHNIQUES  FOR  DETECTING  SUBSURFACE 
CAVITIES  TN  A GRANITE  ENVIRONMENT 

Lewis  S.  Fountain 

Southwest  Research  Institute 


Prepared  for: 

Army  Mobility  Equipment  Research  and 

Development  Center 

Advanced  Research  Projects  Agency 


15  December  1975 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


ffDAOl  82  81 


350Q9  ? 


I EVALUATION  OF  HIGH-RESOLUTION 

EARTH  RESISTIVITY  MEASUREMENT  TECHNIQUES 
FOR  DETECTING  SUBSURFACE  CAVITIES 
IN  A GRANITE  ENVIRONMENT 


by 

Lewis  S.  Fountain 

FINAL  TECHNICAL  REPORT 
Contract  No.  DAAC53-75-C-02 13 
SwRI  Project  14-4250 

Prepared  for 

Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  Va.  22209 

and 

Countermine/Counter  Intrusion  Dept. 

Mine  Detection  Division 
U.S.  Army  Mobility  Equipment 
Research  and  Development  Center 
Fort  Belvoir,  Virginia  22060 

15  December  1975 

Approved  for  public  release;  distribution  unlimited. 

Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

US  P*»pertment  of  Commerce 
Springfield,  VA,  22151 

f 'r[_l 

DEC  12  '915 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  rlWiMi  Data  Fntrrrd) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 


|Z  GOVT  ACCESSION  NO. 


Final  Technical  Report 


« title  imnd  Subiin,)  Evaluation  of  High  Resolution  Earth 
Resistivity  Measurement  Techniques  for  Detecting 
iubsurface  Cavities  in  a Granite  Environment 


7 AUTHOR'D 


Lewis  S.  Fountain 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 RECIPIENT'S  CATALOG  NUMBER 


5 TYPE  OF  REPORT  A PERIOD  COVERFD 

Final  Technical  Report 
12  May  1975-15  Dec.  1975 


* PERFORMING  ORG  REPORT  NUMBER 

Project  14-4250 


■ CONTRACT  OR  GRANT  NUMBERfo 

DAAG5  3-75-C-021  3 


9 PERFORMING  organization  name  and  address 
Southwest  Research  Institute 
P.  O.  Drawer  28510,  8500  Culebra  Road 
an  Antonio,  Texas  78284 


II  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Mine  Detection  Division,  U.  S.  Army  Mobility 
Equipment  Research  and  Development  Center, 


10  PROGRAM  ELEMENT.  PROJECT  TASK 
AREA  A WORK  UNIT  NUMBERS 


1Z  REPORT  DATE 

15  December  1975 


<3  NUMBER  OF  PAGES 


A MONITORING  AGENCY  NAME  A AOORESSCil  di llrrmt  from  Controlling  Other)  IS  SECURITY  CLASS  (of  Hu,  report 

Unclassified 


I5«  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


"6  Distribution  statement  rol  Hum  Report) 


Approved  for  public  release;  Distribution  unlimited. 


*7  OlSTPlBuTlON  5T  ATFMENT  (of  the  mbstrmet  entered  In  Block  20,  II  different  from  Report) 


• 9 KEY  WORDS  f Continue  on  reverse  side  if  necessmry  *nd  identity  by  block  number) 

Earth  resistivity  Borehole  resistivity  measurements 

Pole-dipole  electrode  array  Surface  remote  sensing 

Subsurface  Cavity  detection 
Tunnel  detection  in  granite 

20  ABSTRACT  f Continue  on  reverse  side  II  necessmry  end  identity  by  block  number) 

The  purpose  of  this  work  was  to  experimentally  explore  the  use  of 
earth  resistivity  measurements  using  a pole-dipole  electrode  array  in  the 
detection  of  deep  tunnels  in  a granite  environment.  Measurements  both  on 
the  ground  surface  and  in  drilled  holes  were  evaluated. 

Tests  weie  conducted  at  two  sites  over  existing  mine  adits  or  tunnels, 
one  at  the  Colorado  School  of  Mines  Experimental  Mine  site  and  one  at  a Gold 


EDITION  OF  I NOV  AS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  of  this  PAGE  r»T>rn  Dm.  Entered! 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  This  PAGEf»Fh»n  D«<a  Enfnd) 


O X 3 -m)  cross  section  were  successfully  detected  from  the  surface  at  both 
locations  through  overburden  thicknesses  ranging  from  30  to  HI  feet  (9.  1 to 
24.  7 m).  Other  geological  features  of  the  terrain  were  also  identified  by  the 
measurements.  Successful  detection  was  also  made  from  a dry  borehole 
that  was  lb  feet  (4.  9 m)  horizontally  displaced  from  the  adit  at  the  adit  level. 

Although  near  surface  geological  features  produce  strong  resistivity 
anomalies,  they  can  be  identified.  Small  tunnels  in  granite  can  be  detected 
and  located  to  depths  approaching  100  feet  (30.  5 m)  using  a pole-dipole  surface 
electrode  array.  The  method  is  equally  applicable  for  use  in  boreholes; 
however,  field  conditions  and  drilling  problems  did  not  permit  adequate 
testing  of  the  borehole  resistivity  measurement  technique.  A feature  of 
borehole  measurements  is  that  the  more  1 .mogeneous  environment  below 
the  ground  surface  eliminates  many  of  the  complexities  in  data  interpretation. 
Since  the  sensors  can  be  placed  at  or  near  the  anticipated  tunnel  depth  using 
a borehole,  the  only  maHer  of  concern  is  horizontal  detection  range  from  the 
hole  which  should  at  least  approach  that  of  surface  surveys. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGEOFhan  Dtlt  Enl»r»d) 


SOUTHWEST  RESEARCH  INSTITUTE 
Post  Office  Drawer  28510,  8500  Cuiebra  Road 
San  Antonio,  Texas  78284 


EVALUATION  OF  HIGH-RESOLUTION 
EARTH  RESISTIVITY  MEASUREMENT  TECHNIQUES 
FOR  DETECTING  SUBSURFACE  C(AVIT'ES 
IN  A GRANITE  ENVIRONMENT 


by 

Lewis  S.  Fountain 

FINAL  TECHNICAL  REPORT 
Contract  No.  DA  AC53-75-C-02  I .’I 
SwRI  Proj  oct  I I- 1250 

Prepared  for 

Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  Va  22209 

a nd 

Coinitcrmino/Coiintcr  Intrusion  Dept. 
Mine  Detection  Division 
l .S.  \ r in y Mobility  Equipment 
Research  ami  Development  Center 


5’”  A r-  , ,- 

! n Lo 


nr 

'•  DEC  IS 

Fort  Relvoir.  Virginia  22060  |j  j j-~  - • 


15  December  1975 

Approved  for  public  release;  distribution  unlimited. 


A 


II 


FOREWORD 
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SUMMARY 


The  purpose  of  this  work  was  to  experimentally  explore  the 
use  of  earth  resistivity  measurements  using  a pole-dipole  electrode 
array  in  the  detection  of  deep  tunnels  in  a granite  environment. 

Measurements  both  on  the  ground  surface  and  in  drilled  holes  were 
evaluated. 


Tests  were  conducted  at  two  sites  over  existing  mine  adits 
or  tunnels,  one  at  the  Colorado  School  of  Mines  Experimental  Mine 
site  and  one  at  a Gold  Hill.  Colorado,  site.  The  mine  adits  with  an 
approximate  10  x 10-foot  (3  x 3-m)  cross  section  were  successfully 
detected  from  the  surface  at  both  locations  through  overburden  thick- 
nesses ranging  from  30  to  81  feet  (9.1  to  24.  7 m).  Other  geological 
features  of  the  terrain  were  also  identified  by  the  measurements. 
Successful  detection  was  also  made  from  a dry  borehole  that  was  16 
feet  (4.  9 m)  horizontally  displaced  from  the  adit  at  the  adit  level. 

Although  near  surface  geological  features  produce  strong 
resistivity  anomalies,  they  can  be  identified.  Small  tunnels  in  granite 
can  be  detected  and  located  to  depths  approach'  ig  100  feet  (30.  5 m) 
using  a pole-dipole  surface  electrode  array.  The  method  is  equally 
applicable  for  use  in  boreholes;  however,  field  conditions  and  drilling 
problems  did  not  permit  adequate  testing  of  the  borehole  resistivity 
measurement  technique.  A feature  of  borehole  measurements  is  that 
the  more  homogeneous  environment  below  the  ground  surface  eliminatei 
many  of  the  complexities  in  data  interpretation.  Since  the  sensors 
can  be  placed  at  or  near  the  anticipated  tunnel  depth  using  a borehole, 
the  only  matter  of  concern  is  horizontal  detection  range  from  the  hole’ 
which  should  at  least  approach  that  of  surface  surveys. 
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I. 


INTRODUCTION 


Detecting  and  mapping  of  underground  irregularities  from  the 
earth's  surface  are  of  special  interest  to  military  groups  involved 
in  locating  man-made  tunnels  used  in  guerrilla  warfare.  The  problem 
is  also  of  interest  to  highway  construction  engineers,  archeologists, 
speleologists,  and  mining  and  petroleum  exploration  crews.  As 
important  as  this  problem  is,  no  completely  satisfactory  solution 
has  been  found.  A part  of  the  problem  lies  in  the  complex  structure 
of  the  earth's  thin  aerated  top  soil  which  varies  greatly  in  physical 
‘Jtructure  from  area  to  area  and  even  within  a given  area  making 
exploration  data  difficult  to  interpret.  Geophysicists,  whose  general 
interests  entail  much  greater  penetration  depths,  usually  ignore  the 
first  5-  to  50-foot  (1.5-  to  15.  2-m)  layer  of  soil  by  assigning  a 
correction  factor  to  their  data  to  account  for  errors  involved. 
Consequently,  little  detailed  information  has  been  published  in  the 
literature  concerning  the  small  scale  characteristics  of  the  earth's 
top  layer  of  soil. 

One  characteristic  used  to  classify  the  earth's  substructure  is 
its  apparent  electrical  resistivity  as  me  isured  from  the  surface. 
Measurements  of  earth  resistivity  hav*3  been  used  for  many  years  to 
map  subsurface  stratigraphic  irregularities.  They  can  Drovide  a 
quantitative  measure  of  the  conducting  properties  of  the  subsurface 
and  can  be  used  to  map  the  depth  to  horizons  having  anomalously  high 
or  low  conductivity.  The  method  has  been  used  to  locate  and  map  faults, 
shallow  oil  structures,  gravel  beds,  stratification,  mineralized  faults, 
water  table  depth,  and  to  locate  the  depths  of  transition  from  fresh  to 
salt  water  in  fluid  saturated  formations.  Although  the  reasurement 
method  is  simple,  it  has  not  been  used  extensively  to  explore  the  first 
few  feet  below  the  earth's  surface. 

Southwest  Research  Institute  completed  a series  of  field  tests 
in  June  1975  to  evaluate  a high- resolution  oole-dipole  earth  resistivitv 
survey  technique  for  detecting  and  mapping  sinkhole  cavities  in 
connection  with  highway  routing  and  maintenance  problems1.  Analysis 
of  the  test  results  obtained  with  this  surface-operated  geophysical 

L.  S.  Fountain,  F.  X.  Herzig,  and  T.  E.  Owen,  "Detection  of  Subsurface 
Cavities  by  Surface  Remote  Sensing  Techniques,"  Report  No.  FHWA-RD- 
75-80,  SwRI  Final  Technical  Report,  Contract  No.  DOT-FH- 1 1 -8496, 
Southwest  Research  Institute  Project  14-4250,  June  1975. 
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measurement  method  showed  excellent  success  in  locating  small  under- 
ground solution  cavities  in  Florida  limestone  and  in  locating  subsurface 
cavities  and  other  geologic  anomalies  along  highway  rights  - or- way  in 
Alabama  and  Florida.  A graphical  analysis  was  used  that  allowed  the 
position  of  the  detected  voids  to  be  determined. 

This  same  method  had  been  successfully  used  to  locate  solution 
cavities  in  limestone  environments  by  Bates2  and  Bristow3. 

A problem  of  concern  to  the  United  States  Army  is  the  detection 
of  cavities  and  tunnels  in  granite  at  depths  of  from  a few  feet  (about 
1 metre)  to  greater  than  150  feet  (45.  7 m)  below  the  ground  surface 

to'loVlTfT/f  tlleSVargetS  are  in  the  ranSe  6 x 6 feet  (1.  8 x 1.  8 m) 

, *:e  ( x 3 m).  These  depths  and  target  dimensions  were  in 

e range  o detection  capabilities  of  the  earth  resistivity  method,  but 
such  measurements  had  not  been  evaluated  in  a granite  environment. 

It  was  the  purpose  of  this  program  to  evaluate  earth  resistivity 
in  granite ^ " * POSSlble  solutio"  to  the  detection  of  deep  tunnel. 

The  scope  of  the  program  covered  the  buildup  of  simple  instru- 

To!h  atlh"  f°;  feld  data  Collection;  collection  of  earth  resistivity  data 

and  sin  b°reh°le^.and  on  the  earth's  surface  at  two  test  sites  in  Colorado; 
and  simple  graphical  analysis  of  the  collected  da.^a 


E.R  Bates  "Detection  of  Subsurface  Cavities,"  Miscellaneous 

' a Pe^  ~ 762538,  U.  S.  Army  Engineer  Waterways  Exper- 

iment Station,  Vicksburg,  Mississippi,  June  1 973. 

LrCF5nrif°W’  ''.A  NeW  GraPhical  Resistivity  Technique  for  Detecting 
Decembe^  1%6  ' " ^ b Part  4.  204-227, 
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II. 


EARTH  RESISTIVITY  SEARCH  METHOD  USING  THE 
POLE-DIPOLE  ELECTRODE  ARRAY 


P 


A.  General  Background 

A variety  of  electrical  resistivity  geophysical  explo- 
ration methods  have  been  explored  as  possible  approaches  to  sub- 
surface cavity  detection.  The  basic  resistivity  measurement  and 
analysis  concept  underlying  these  various  attempts  has  been  essentially 
the  same  in  each  case:  that  is,  by  establishing  an  otherwise  predictable 
electrical  current  distribution  within  a relatively  Large  volume  of 
homogeneous  earth  material,  any  observed  perturbations  in  the  current 
distributions  measured  as  potentials  or  electr.'c  fields  at  the  ground 
surface  can  be  interpreted  in  terms  of  possible  subsurface  structural 
or  earth  material  resistivity  anomalies.  The  degree  of  perturbation 
in  the  current  distribution  is  dependent  upon  the  resistivity  contrast 
between  the  anomalous  subsurface  structures  and  the  surrounding  earth 
material  and,  equally  important,  the  detectability  of  such  perturbations 
is  also  dependent  upon  the  size  and  shape  of  the  anomaly  and  its  orienta- 
tion relative  to  the  current  flow. 

The  main  differences  among  the  various  electrical 
resistivity  geophysical  profiling  methods  are  largely  in  the  electrode 
array  patterns  used  to  establish  the  subsurface  current  distributions 
and  in  measuring  the  potential  differences  at  the  ground  surface.  The 
manner  in  which  the  electrodes  are  moved  or  scanned  over  the  area 
being  surveyed  also  differs  with  the  different  electrode  arrays  as  do 
the  methods  of  resistivity  data  analysis  and  interpretation. 

Detailed  descriptions  and  discussions  of  various  electrical 
resistivity  methods  of  subsurface  cavity  detection  are  presented  in 
Appendix  A. 

The  pole-dipole  array  has  had  less  previous  use  and  has 
shown  the  greatest  potential  for  detecting  underground  cavities  and 
predicting  their  depths  and  locations.  This  method  will  be  discussed 
in  detail  below. 

B.  Pole-Dipole  Earth  Resistivity  Electrode  Array 

1.  Theory  and  Method 

The  pole-dipole  electrical  resistivity  survey 


i 
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method  is  based  on  a tour-electrode  array  configuration  in  which 
the  current  sink  electroce  is  located  at  effective  infinity,  and  the 
potential  electrodes  are  separated  from  one  another  by  a fixed 
minimum  distance  proportional  to  the  desired  resolving  power  of  the 
system.  The  potential  electrode  pair  is  located  at  various  positions 
a ong  a selected  line  near  the  current  source  electrode  as  the  means 
of  vertically  sounding  the  subsurface  below  that  electrode.  The 
current  source  electrode  is  then  moved  ahead  along  the  line  at  suitable 
incremental  distances  to  provide  Horizontal  profile  scanning  The 
pole-dipole  array  is  illustrated  in  Figure  1.  In  order  that  the 


Hernia  phe  rical 
Equipotential  Surfaces 
And  Resistivity  Resolution 
Volume 


Current 

Source 

Electrode 


FIGURE  1 . POLE- DIPOLE  EARTH  RESISTIVITY 
ELECTRODE  ARRAY 


equi potential  surfaces  be  hemispherical  ard  concentric  about  the 
source  electrode,  the  sink  electrode  must  be  located  at  effective  infinity 
w ich  will  generally  be  about  5 to  10  times  the  largest  value  of  detection 
penetration  depth  of  interest  in  the  survey.  In  a typical  survey  the 

potential  electrode  pair  is  spaced  at  a fixed  separation  of  20  feet  (6  1 rr 

nd  moved  away  from  the  source  electrode  at  10-foot  (3-m)  intervals  to 
provide  overlapping  and  redundant  resistivity  data. 
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. . , A modified  form  of  this  array  has  been  used 

:;Z“e  r0!)‘a,biUl!'  and  accurac>'  of  detecting  email  subsurface 
to  edifications  pertain  largely  to  the  collection  of  over- 
pping  field  data  using  closer  potential  electrode  spacings  These 
overlapping  measurements  are  elaborate  enough  to  provide  up  to  four 

for  rea  em  t , theleCtr0dePair  readin«s  at  each  measurement  station 
or  refinement  of  the  resistivity  interpretation  analysis.  This  modi- 

Thta  ,°.\h e"va‘]S  meaaurements  on  both  sides  of  the  source  electrode. 
This  is  the  method  used  by  SwRI  in  the  work  reported  in  Reference  1 and 

tannin.  COVered  *7  ^ F°r  clarity  °‘  ^P^affon.  the  over- 

lapping  survey  procedure  is  described  as  follows  for  a typical  100-foot 

(:>J.  b-m)  penetration  depth  survey: 


a*  Place  the  current  source  electrode.  Ci 
as  shown  in  Figure  1.  at  the  first  traverse  station; 

. . b*  Place  the  current  sink  electrode,  C? , at 

a minimum  ci^tauce  of  500  feet  ( i 52.  4 m).  (5  r2),  behindC.  on  the  pre- 
s ablished  traverse  line  having  already  decided  that  the  maximum 
potential  electrode  scan  distance  from  Cl  will  be  100  feet  (30.  5 m); 

c-  Place  potential  electrodes,  PI  and  P2  at 

pre-established  station  markers  at  10  feet  (3  m)  and  20  feet  (6.  1 m) 

respectively,  on  the  sink  electrode  side  of  Cl,  and  obtain  resistance 
resding  j 


7 ^ doove  wicn  potential  elect 

30  f'  ,q\  at  Pre-es.ablished  station  markers  at  20  feet  (6.  1 m)  and 

and  “as  m,  **  * 9°  f“‘  <”•  * “ 


e‘  Next,  place  potential  electrodes  at  pre- 
es  a isned  10-foot  (3-m)  interval  station  markers  on  the  opposit^  side 

° inTfQQ,a/  n resi8tance  ladings  at  each  station  pair  out 

for  th  f t f"m)  limit'  ™S  comPletes  the  survey  procedure 

or  the  first  current  station  position  for  Cl ; 


of  40  feet 
(30.  5-m) 


f-  M°ve  Ci  up  the  traverse  line  a distance 

m),and  obtain  resistance  readings  over  the  100-foot 
scan  zones  on  each  side  of  this  current  station; 
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Repeat  step  f above  until  complete  traverse 

The  above  survey  procedure  gives  four  overlapping  resistance  readings 
for  each  10-foot  (3-m)  spaced  potential  electrode  pair  station  to  a depth 
of  more  than  ^0  feet  (IS,  2 m)  as  the  survey  proceeds.  This  procedure 
with  both  5-  and  10-foot  (l.  5-  and  3-m)  potential  electrode  spacing  was 
used  in  all  of  the  pole-dipole  earth  resistivity  tests  performed  cn  this 
program.  The  bi-directional  potential  scanning  doubles  the  amount  of 
data  per  current  station  over  that  obtained  by  earlier  procedures.  The 
four-level  data  overlap  aids  in  resolving  the  locations  and  resistivity 
contrasts  of  the  various  high-  and  low- resistivity  anomalies  which  may 
be  encountered  along  the  traverse.  This  is  accomplished  by  the  multiple 
positional  aspect  observations  of  possible  anomalies  as  detected  from 
several  survey  stations;  whereas,  in  contrast,  for  less  redundant  obser- 
vations a high- resistivity  anomaly  and  a low- resistivity  anomaly 
occurring  in  the  same  hemispherical  shell  field  of  view  will  tend  to 
nullify  one  another  and  be  undetected. 


g- 

line  is  surveyed. 


A straightforward  graphical  analysis 
method  was  devised  capable  of  utilizing  and  displaying  all  of  the  field 
data  as  a means  for  locating  the  experimental  best-fit  positions  and  depths 
of  detected  subsurface  cavities.  The  success  of  this  analysis  approach 
is  largely  achieved  through  the  spatial  redundancy  of  the  field  data  with 
the  result  that  target  ambiguities  and  false  interpretations  are  minimized, 
and  improved  cavity  size  and  shape  indications  are  derived. 


Graphical  Data  Analysis 


The  graphical  analysis  of  the  pole-dipole  resistivity 
data  can  best  be  illustrated  by  first  examining  a sample  of  field  data.  The 
basic  field  measurements  as  performed  on  this  program  were  recorded  on 
specially  prepared  data  forms.  Figure  2 is  a sample  data  sheet  showing 
recorded  instrument  readings,  calculated  resistivity  values,  and  a graph 
of  a derived  resistivity  profile.  The  first  column  on  the  data  sheet  is  the 
potential- pair  electrode  distance  from  the  current  source  electrode;  the 
second  column  lists  the  distances  from  current  electrode  to  midpoint 
between  potential  electrodes;  the  third  column  is  the  resistivity  instru- 
ment reading  in  ohms  (or  in  volts  if  AC  instrumentation  if  used);  the 
fourth  column  is  the  geometrical  factor  required  to  calculate  apparent 
resistivity  for  the  distances  listed  in  the  first  column  using  the  resistances 
recorded  in  the  third  column  or  calculated  from  the  recorded  voltage  and 
eart  i current  for  AC  instrumentation;  and  the  fifth  column  is  the  calculated 
value  of  apparent  resistivity  in  ohm  - centimete  rs  x 103  (S2  • m x 10). 
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The  pole-dipole  geometrical  factor  in  column 
four  is  calculated  from  the  relationship, 

Z-nr\rz 

K = (1) 

rz  - n 

where  r^  and  r 2 are  the  distances  of  the  potential  electrodes,  Pj  and  P2. 
from  the  current  source  electrode,  Cj,  as  illustrated  in  Figure  1.  The 
basic  pricniples  on  which  the  geometrical  factor,  P,  is  derived  are 
discussed  in  Appendix  A and  applied  to  several  earth  resistivity  electrode 
array  configurationa 


Apparent  resistivity  is  calculated  from  the  equation 
developed  in  Appendix  A 

p = KR  (2) 

where  R is  the  measured  resistance. 

A graph  of  the  apparent  resistivity  profile  on  one 
side  of  the  current  source  electrode  is  plotted  on  the  data  sheet  as  shown 
in  Figure  2,  Those  points  on  the  profile  that  indicate  high  or  low  resis- 
tivity perturbations  away  from  the  average  profile  trend  are  next 
identified  and  marked  for  transfer  to  a scaled  drawing  used  to  graphically 
locate  the  anomalous  underground  resistivity  structures.  Determinations 
of  the  average  profile  trends  are  made  by  visual  inspection.  A simple 
and  liberal  interpretation  of  high  and  low  resisti  dty  perturbations  is 
permissible  since  the  redundancy  of  the  data  anc  a required  multiplicity 
of  pe  turbations  associated  with  each  underground  anomaly  will  delete 
the  improper  interpretations. 

An  example  of  the  scaled  drawings  used  in  the 
graphical  analysis  is  set  up  as  shown  in  Figure  3.  Thu  distance  along 
the  ground  surface  representing  the  survey  traverse  is  marked,  and  the 
consecutive  positions  of  the  current  probe  are  shown  by  the  arrows. 

The  perturbations  interpreted  from  the  forward  and  reverse  resistivity 
profiles  for  each  current  electrode  station  are  denoted  on  the  bracketed 
lines  drawn  above  the  ground  surface  line.  The  high  resistivity  anomalies 
in  this  example  are  labeled  as  A,  B,  C,  and  D as  taken  from  the  related 
data  graphs  shown  at  the  top  of  the  illustration.  The  data  in  this 
example  are  actual  field  data  that  show  detection  of  a 10  x 10-foot 
(3  x 3-m)  mine  adit  30  feet  (9.  1 m)  below  the  surface.  With  a compass 
centered  at  each  current  location  on  the  ground  surface  line,  arcs  are 
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drawn  at  distances  representing  the  bounds  of  the  high  resistivity 
anomaly  with  respect  to  the  current  probe  position.  The  pairs  of  circular 
arcs  are  labeled  A,  B,  C,  and  D corresponding  to  the  interpreted 
< nomalies.  The  space  where  the  four  sets  of  arcs  intersect  is  the 
graphically  derived  location  of  the  underground  structure  responsible 
for  the  high  resistivity  perturbations.  Indeed,  a useful  arbitrary  guide 
for  taking  advantage  of  the  redundancy  of  the  overlapping  field  measure- 
ments is  to  require  a minimum  of  three  arcs  to  intersect  at  a common 
anomaly  location  before  it  is  interpreted  with  any  confidence  as  a 
probable  underground  cavity.  Moreover,  as  illustrated  in  Figure  3 
the  arcs  are  drawn  only  in  the  90-degree  sectors  corresponding  to  either 
the  forward  or  reverse  profiles  containing  the  perturbations  being  used. 

The  reason  for  this  is  that  the  distortions  of  the  equipotential  lines 
(represented  in  a first  order  manner  by  the  circular  arcs)  are  very  weak 
if  the  perturbing  anomaly  is  located  in  the  opposite  bO-degree  sector. 

It  should  be  recognized  that  the  pa^rs  of  arcs  drawn 
for  each  resistivity  perturbation  not  only  describe  the  resistivity  pertur- 
bations in  the  vertical  plane  along  the  ground  traverse  line  but  applies, 
to  some  extent,  to  a three-dimensional  spherical  shell  segment  extending 
laterally  on  each  side  of  the  traverse  line.  It  has  been  speculated  by  both 
Bristow  and  Bates  that  this  apparent  lateral  *ield  of  view  away  from  the 
traverse  line  is  contained  within  an  angle  of  auout  + Zb  degrees  on  each 
side  of  the  traverse  line  relative  to  the  source  current  electrode  location. 

Detailed  interpretation  of  tl  e subsurface  arc  inter- 
sections is  complex,  and  a knowledge  of  the  geology  of  the  area  is  almost 
a necessity.  In  many  cases  the  high  resistivity  anomalies  are  caused  by 
bedrock-soil  interfaces,  especially  in  locations  where  the  bedrock  elevations 
are  very  irregular.  A low  resistivity  indication  is  usually  a soil-  or 
solution-fsoil  with  a very  high  water  content)  filled  cavity  surrounded  by 
otherwise  higher  resistivity  structures  such  as  dolomite  or  granite.  A 
few  test  borings  in  the  test  area  can  be  a great  help  in  interpreting  the 
detailed  implications  of  the  anomalies. 

Errors  in  data  interpretation  can  be  caused  by 
underground  wires,  metal  fences  with  metal  posts  on  the  surface,  buried 
pipes,  abrupt  cliffs,  ponds,  streams,  etc.  in  the  area  surveyed.  The 
effects  of  these  various  disturbing  features  is  to  distort  the  ee  rth  current 
patterns  from  the  assumed  spherically  divergent  paths  postulated  for  a 
flat  homogeneous  ground.  Interference  features  of  this  type  should  be 
avoided  where  possible  or  approached  with  the  electrode  array  oriented 
at  right  angles  with  re  jpect  to  the  obstacle. 
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3. 


Application  in  Boreholes 


a.  Method 


In  applying  the  method  in  a single  vertical 
borehole,  a resistivity  measurement  cable  string  consisting  of  a current 
source  electrode  and  two  potential  electrodes  is  located  in  the  borehole; 
the  current  sink  electrode  is  located  at  effective  infinity  on  the  surface 
as  was  done  in  the  surface  survey  measurements.  Measurements  are 
then  made  by  moving  the  fixed- space  pair  of  potential  electrodes  in 
chosen  increments  away  from  the  current  source  electrode.  Several 
sets  ot  measurements  can  be  made  in  a single  borehole.  The  current 
source  electrode  can  first  be  placed  in  the  bottom  of  the  hole  and 
resistivity  measurements  made  upward  from  the  bottom  of  the  hole  toward 
the  top.  The  current  source  electrode  is  then  placed  at  another  location 
up  the  hole  and  resistivity  measurements  made  around  this  new  current 
source  location.  This  technique  is  repeated  up  the  hole  until  the 
measurements  are  well  above  the  anticipated  depth  of  the  target  anomaly. 

. he  sets  of  data  can  then  be  used  with  the  same  graphical  method  as  is 
used  with  surface  surveys  to  obtain  a depth  below  surface  bearing  on  a 
resistivity  anomaly  such  as  might  be  caused  by  a subsurface  void. 

When  the  method  is  applied  below  the 
ground  surface  in  a homogeneous  earth,  the  current  injected  at  the 
source  electrode  will  diverge  with  radial  symmetry  into  the  earth 
medium  giving  rise  to  concentric  spherical  equipotential  surfaces 
equivalent  to  the  hemispherical  surfaces  described  earlier.  Resistivity 
anomalies  occurring  within  any  two  closely  spaced  spherical  equipotential 
surfaces  (a  potential  difference  shell)  will  be  observed  although  the 
azimuthal  orientation  of  such  anomalies  will  not  be  discernable  until 
supplemental  measurements  are  made  in  other  nearby  boreholes.  From 
a single  borehole,  however,  anomalous  geological  features  can  be 
detected,  resolved  in  size,  and  located  accurately  in  depth  through  the 
process  of  making  overlapping  and  redundant  resistivity  measurements 
along  the  borehole.  This  concept  is  illustrated  in  Figure  4 where  the 
anomalous  geologic  target  occupies  only  a localized  volume.  Localization 
of  the  anomaly  in  depth  is  evident  from  the  intersection  of  the  two 
spherical  potential  shells  shown;  however,  the  intersection  zone  of  these 
shells  is  toroidal  in  shape  and  oriented  symmetrically  about  the  borehole 
axis.  Thus,  if  the  aximuthal  position  of  the  target  relative  to  the  borehole 
axis  is  to  be  determined  without  ambiguity,  it  must  be  done  by  the  inter- 
section of  at  least  three  ether  similar  toroidal  detection  zones  responding 
to  the  same  target  as  derived  from  resistivity  measurements  made  in 
separate  adjacent  boreholes. 


1 l 


Surface 


FIGURE  4.  SIMPLIFIED  ILLUSTRATION  OF  SINGLE-BOREHOLE 
HIGH  RESOLUTION  RESISTIVITY  MEASUREMENT  METHOD 


The  geometrical  factor  used  in  the 
calculation  of  apparent  resistivity  from  measurements  made  in  a 
borehole  changes  when  the  current  source  is  below  the  surface  to 

47Tri  r2 

kb  — u) 

This  change  is  because  the  potential  difference  shells  around  the 
current  source  electrode  are  now  spherical  in  shape. 

b.  Detection  Parameters 


Many  factors  affect  the  underground 
detection  effectiveness  of  any  earth  resistivity  search  method.  Some 
of  the  more  important  factors  are: 


(1) 

(2) 

(3) 

line  (distance  from  borehole  or,  for 
surface) 


Target  shape 
Target  size 

Target  distance  from  traverse 
surface  survey,  depth  below 


source 


(4)  Target  distance  from  current 


target  and  background 


(5)  Potential  probe  spacing 

(6)  Resistivity  contrast  between 


A study  done  by  SwRI  for  the  U.  S.  Bureau 
of  Mines  on  another  program  has  evaluated  some  of  these  parameters, 
as  they  relate  to  resistivity  measurements  from  a borehole.  A general 
summary  of  the  results  will  be  given  here  to  provide  some  insight  into 
detection  possibilities. 


Target  shape  is  a critical  factor  as  is 
its  size.  The  distance  to  which  a spherically  shaped  target  may  be 
defected  is  directly  proportional  to  its  diameter  and  inversely  pro- 
portional to  the  cube  root  of  the  certainty  with  which  the  background 
signal  can  be  predicted.  If  the  geological  background  is  sufficiently 
uniform  that  the  background  resistivity  signal  can  be  predicted  to 
within  1%,  then  signal  perturbations  of  2%  should  Le  recognizable. 
Under  such  conditions  a spherically- shaped  target  of  good  resistivity 
contrast  could  be  detected  at  a distance  of  almost  3 radii  from  the 
borehole  to  the  centroid  of  the  target  (or  about  one  target  diameter 
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distance  from  borehole  to  nearest  surface  of  the  target).  A cylindrically 
shaped  target  having  the  same  strength  perturbation  signal  as  the  sphere 
would  he  detectable  at  a greater  distance  than  the  sphere.  For  example, 
if  a 10-foot  (3-m)  diameter  spherical  void  centered  15  feet  (4.  5 m)  from 
the  borehole  produces  a given  signal  above  noise,  then  a long  cylindrical 
cavity  with  the  same  diameter  would  produce  the  same  signal  at  approxi- 
mately 53  feet  (16.  2 m).  This  relationship  is  expressed  as 


Distance  to  Cvlind^r  - K /Distance  to  Sphere  ,1/2  . 

n j- ) Distance  to  (4) 

Sphere  Radius  17 

Sphere 

where  K is  a function  of  the  ratio  of  target  and  background  conductivity. 
In  the  case  of  an  air-filled  void,  K = 2. 


. . . Potential  electrode  spacing  used  in  borehole 

resistivity  measurements  plays  an  important  part  in  resolving  targets 
Numerical  modeling  has  shown  that  for  electrode  spacing  less  than  half 
the  arget  sphere  ^diu8,  the  detection  capability  is  essentially  independenl 
e ec  rode  spacing.  For  larger  spacing  the  sensitivity  falls  off  linearly 
as  pro  e spacing  increases.  As  an  example,  the  perturbation  component 
of  a processed  signal  for  an  electrode  spacing  equal  to  two  target  sphere 
radii  is  about  one  half  its  largest  value.  H 


it  is  intuitive  that  an  air-filled  void  within 
a omogeneous  volume  of  conducting  earth  material  would  impose  a 
greater  resistance  to  electric  currents  flowing  in  the  vicinity  of  the 
void  than  in  other  parts  of  the  material  volume.  Conversely,  a water- 

t°hem  H r fV°ldf  TU  °ffer  a l°Wer  resistance  to  electric  current  since 
e conductivity  of  the  filling  material  will  generally  be  greater  than  that 

of  the  surrounding  material.  The  degree  of  perturbation  in  the  current 

distribution  is.  therefore,  obviously  dependent  upon  the  resistivity 

contrast  between  the  anomalous  subsurface  structures  am',  ^he  surrounding 

from  T'T  8,  • ^ analysis  showed  that  the  ability  to  detect  any  target 

CTar,etr'  H I " t0  **  COnd^ivity  contrast  between 

the  target  and  the  surrounding  materials. 


C. 


Instrumentation 


DC  Instrument 


Earth  resistivity  measurements 
both  AC  and  DC  type  instruments. 


were  made  using 
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The  DC  measurements  were  made  using  a Keck 
Model  1C-69  instrument  aj  shown  set  up  in  the  field  along  with  acces- 
sories in  Figure  5.  Figure  6 is  a photograph  showing  two  metal  current 
electrodes  and  two  porous-pot  potential  electrodes  used  with  the  instru- 
ment. 

The  commercial  instrument  is  a DC  system 
obta  ning  power  for  earth  current  from  dry  cell  batteries.  Four  45V 
batteries  are  housed  in  the  instrument  case,  and  there  is  a receptacle 
on  the  case  where  additional  batteries  can  be  added  \s  required.  The 
battery  case  shown  on  tie  ground  to  the  right  of  the  instrument  in  Figure 
^-c”tains  10  additional  NEDA  205,  45V  batteries.  Switches  on  the 
instrument  panel  and  on  the  battery  pack  allow  individual  batteries  to 
be  switched  into  the  circuit  in  series  as  required.  Resistance  can  be 
measured  with  the  instrument  over  the  range  0.  001-1000  ohms,  and 
the  dial  can  be  read  to  one  part  per  thousand. 

The  Keck  resistivity  instrument  can  also  be  used 
to  measure  self  potential  that  often  results  from  the  electrochemical 
action  between  soil  solutions  and  ore  bodies.  In  nearly  all  cases  when 
two  electrodes  are  placed  on  the  ground  surface  a DC  potential  can  be 
observed  between  them.  If  porous-pot  electrodes  are  used,  this 
potential  will  be  related  to  natural  earth  currents.  The  instrument 
co’  tains  a circuit  capable  of  nulling  out  this  natural  earth  potential 
so  that  it  does  not  affect  the  resistivity  measurements.  The  instru- 
ment dial  used  to  adjust  the  indicator  circuit  for  a null  is  calibrated 
to  read  in  millivolts  and  can  be  read  when  a null  is  obtained.  This 
self  potential  data  can  be  of  value  in  aiding  interpretation  of  earth 
resistivity  measurements  collected  over  solution-filled  cavities. 

The  instrument  has  a "forward- reverse"  switch  on  the  panel  to 
change  the  direction  of  the  current  flow  into  the  earth.  For  precise 
resistance  determinations,  field  measurements  were  made  using  both 
directions  of  current  flow,  and  the  two  measurements  averaged. 

Another  reversing  switch  was  added  as  a modi- 
fication to  the  instrument.  This  was  done  to  allow  faster  setup  of 
electrodes  in  special  electrode  arrays.  The  potential  electrodes  must 
be  connected  to  the  instrument  in  the  proper  polarity  since  the  unit  is  a 
DC  instrument,  and  DC  nulling  potentials  are  used.  Once  in  operation 
if  the  two  potential  electrodes  are  reversed,  as  might  occur  if  alternate 
electrodes  are  moved  ahead  along  the  survey  traverse,  the  additional 
reversing  switch  allows  the  potential  electrodes  to  be  reversed  inde- 
pendently of  the  current  electrodes. 
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FIGURES.  EARTH  RESISTIVITY  INSTRUMENT 
AND  ACCESSORIES  SET  UP  FOR  OPERATION  IN  THE  FIELD 
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FIGURE  6.  METAL  CURRENT  ELECTRODES 
AND  POROUS-POT  POTENTIAL  ELECTRODES  USED 
WITH  THE  EARTH  RESISTIVITY  INSTRUMENT 


Porous-pot  electrodes  were  used  as  potential 
electrodes,  and  metal  electrodes  were  used  as  current  electrodes. 
Porous-pot  electrodes  are  non-polarizing  electrodes  used  to  eliminate 
problems  and  errors  caused  by  galvanic  action  between  metal  electrodes 
and  the  soil. 


Electrical  contact  with  the  earth  using  a porous-pot 
is  made  through  the  moist  bottom  of  the  pot.  The  body  of  the  pot  is 
glazed  porcelain  except  for  the  pottom.  Saturated  copper  sulphate 
solution  inside  the  pot  permeates  the  porous  ceramic  bottom  to  form  the 
electrical  ground  surface  contact.  A pure  copper  electrode  immersed 
in  the  solution  and  passing  through  the  cover  on  the  pot  permits  con- 
nection to  the  instrument. 

Pusher  style  metal  electrodes  as  shown  in  Figure  6 
were  used  as  current  electrodes  except  in  rocky  areas  where  electrodes 
must  be  driven  in.  The  current  electrode  assembly  is  4 feet  (1.  2 m) 
long  including  a 1-foot  (0.  3-m)  long,  1 /2-inch  (13-mm)  diameter  copper 
clad  steel  electrode  at  one  end  attached  to  a rod  having  a horizontal 
extension  that  allows  an  operator  to  push  it  into  the  soil  using  foot 
pressure.  No  ground  contact  problems  were  encountered  using  the 
metal  current  electrodes. 

The  cable  reel  set  shown  in  Figure  5 provided  a 
convenient  means  for  handling  the  wire  leads  used  with  the  current 
and  potential  electrodes.  This  reel  assembly  contains  four  reels  of 
vinyl  insulated  18  gauge  wire,  each  containing  1000  feet  (305  m)  of 
wire.  The  cable  reels  are  supported  in  bronze  bearings  and  have  hand 
cranks  for  rapid  rewinding  of  the  wire.  Connections  to  the  instrument 
are  made  through  copper  disks  mounted  on  each  reel  and  through  heavy 
carbon  brushes  contacting  the  copper  disks.  The  brushes  connect  to 
a quick  disconnect  connector  on  the  outside  of  the  case.  Cable  reels 
are  insulated  from  the  case.  A short  jumper  cable  connects  the  reel 
case  to  the  instrument. 

2,  AC  Instrument 


Alternating  current  resistivity  measurements 
were  also  used  for  two  reasons:  (1)  surveys  can  be  made  much  faster 
because  there  is  no  requirement  to  balance  out  DC  earth  materials, 
and  (2)  the  borehole  measurements  necessarily  used  metallic  potential 
electrodes  which  would  give  rise  to  galvanic  potentials  in  a DC  system. 
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The  AC  instrumentation  was  laboratory  constructed 
and  comprised  a constant  current  power  source  and  a high  input  impedance 
autoranging  voltmeter  with  digital  readout. 


The  AC  current  source  operated  at  a frequency  of 
97  Hz  to  prevent  signal  interference  from  power  line  frequencies  and 
harmonics.  The  constant  current  capability  was  necessary  to  prevent 
electrode  contact  resistance  variations  from  causing  errors. 

The  input  signals  to  the  voltmeter  were  amplified 
using  a differential  amplifier  arrangement  and  narrowband  filtered  to 
eliminate  noise  before  it  was  connected  to  the  voltmeter  circuit.  The 
input  impedance  was  in  the  range  of  25  megohms,  thus,  keeping  errors 
caused  by  potential  electrode  contact  resistance  variations  to  a minimum 


A great  advantage  in  using  AC  instrumentation  for 
surface  resistivity  measurements  is  faster  measurements.  With  the  DC 
instrument,  time  is  required  to  first  adjust  a self  potential  dial  for  a 
null  reading  and,  then,  to  adjust  another  dial  for  a null  reading  to  obtain 
the  desired  resistance  value.  When  AC  power  is  used,  the  polarizing 
effects  around  electrodes  are  eliminated  by  the  current  flow  reversals. 
Therefore,  since  most  unwanted  ground  currents  are  DC  no  compen- 
sation or  nulling  adjustment  is  required.  Using  an  autoranging  voltmeter 
to  read  potentials,  from  which  resistance  is  finally  calculated,  the  reading 
for  the  measurement  can  be  made  immediately  after  the  electrodes  are 
placed  in  contact  with  the  earth. 


Electrodes 

all  metal.  Stainless  steel,  steel, 
all  used  successfully. 


used  with  AC  instrumentation  were 
copper  clad  steel,  and  aluminum  were 


Special  sPring  potential  electrodes  were  fabricated 
for  use  in  the  boreholes.  Both  berylium  copper  and  stainless  steel  were 
used  to  construct  these  electrodes  as  shown  in  Figure  7.  The  electrodes 
were  assembled  on  a special  section  of  one-inch  (25.  4-mm)  PVC  pipe. 
Additional  10-foot  (3-m)  sections  were  fabricated  with  special  end 
connectors  so  that  the  potential  measuring  pair  of  electrodes  could  be 
lowered  by  hand  to  any  depth  in  the  borehole. 

A short  cylindrical  aluminum  electrode  was 
constructed  to  be  used  as  a current  electrode  that  could  be  lowered 
to  the  bottom  of  a borehole  to  make  electrical  contact  through  a small 
amount  of  water  poured  into  the  hole.  The  electrode  is  approximately 
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FIGURE  7.  PHOTOGRAPH  OF  A BOREHOLE  RESISTIVITY 


8 inches  (20  3 mm)  long,  1-1/4  inches  (32  mm)  in  diameter,  and  weighs 
4-1  2 pounds  (2  kg).  It  was  made  by  filling  a section  of  aluminum  tubing 
with  lead  and  attaching  an  electrical  terminal  and  an  eye  connection  for 
mechanical  attachment  for  raising  and  lowering  it  in  a hole. 

A special  waterproof  cable  was  made  up  for  use 
with  the  potential  electrodes.  It  is  a Belden  8728  cable  with  a poly- 
urethane jacket.  The  cable  is  20  feet(6.  1 m)  long  terminated  at  the 
electrode  end  by  an  encapsulated  hign-input  impedance  preamplifier 
having  waterproof  leads  that  attach  to  the  potential  electrodes.  The 
other  end  of  the  cable  connects  to  a filter  amplifier  and  voltmeter  at 

the  ground  surface.  The  cable  was  made  waterproof  for  anticipated 
future  wet-hole  use. 


III. 


FIELD  TESTS 


A.  Idaho  Springs , Colorado,  Site 
1 . Site  Description 

The  test  site  at  Idaho  Springs,  Colorado,  was  at 
the  Colorado  School  of  Mines  Experimental  mine.  Located  at  an 
elevation  of  about  8000  feet  (2.  4 km)  the  mine  adit  extends  into  the 
mountain  a distance  of  about  5 miles  (8  km).  Its  cross  section  is 
aoproximately  10  x 10  feet  (3x3  m)  cut  into  almost  solid  granite. 

The  mountain  slopes  upward  from  the  adit  entrance  at  a 30°  angle 
making  it  possible  to  selectany  desired  overburden  thickness  when 
field  tests  are  made.  The  very  thin  soil  is  a sandy  composite  of 
mostly  weathered  granite.  A good  portion  of  the  surface  area  above 
the  adit  is  solid  rock.  Figure  8 is  a view  of  the  test  site  showing  the 
entrance  to  the  adit  and  the  very  rough  steep  slope  where  resistivity 
measurements  were  made.  Figure  9 is  a photograph  showing  the 
sloping  area  and  the  type  of  surface  encountered  during  the  tests. 

Figure  10  is  another  photograph  taken  in  the  area  showing  the  rock 
that  comes  up  to  the  surface.  The  hole  is  a shaft  going  into  the  mine. 

2.  Site  Preparation 

The  entire  area  over  the  mine  adit  that  was  to 
be  used  in  the  tests  was  surveyed,  and  contour  maps  were  made. 

This  was  done  to  aid  in  surface  resistivity  data  interpretation  and  to 
locate  a borehole  in  the  road  that  was  to  be  used  for  borehole  resis- 
tivity measurements.  Figure  1 I is  a contour  map  of  the  entire  area 
used,  with  contour  intervals  of  20  feet  (6.  I m).  From  this  map  traverse 
lines  were  selected  to  allow  resistivity  surveys  to  be  made  across  the 
adit  where  it  was  50  and  73  feet  (15.  2 and  22.  2 m)  below  the  surface. 
Figure  12  is  a contour  map  of  the  surface  from  the  adit  level  tc  the  service 
road  that  crossed  over  it  at  an  elevation  of  30  feet  '9.  I m).  The  Sections 
shown  on  this  map  indicate  sections  through  a borehole  that  will  be 
discussed  later. 


Layout  plans  included  two  boreholes  to  be  put 
into  the  service  road  that  would  miss  the  adit  by  distances  of  16  and  40 
feet  (4.9  and  12.  2 m).  Drilling  difficulties  prevented  completion  of  one 
of  tne  holes.  Only  the  one  that  had  a mis  s - distance  of  16  feet  (4.  9 m) 
was  completed  and  used  in  the  tests.  It  was  drilled  to  a depth  of  60 
feet  (18.  3 m)  which  was  about  20  feet  (6.  1 m)  deeper  than  the  floor  of 
the  adit.  Heavy  rains  after  the  hole  was  completed  washed  about  10  feet 
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FIGURE  8.  VIEW  OF  THE  COLOR  ADO  SCHOOL  OF  MINES 
TEST  SITE  SHOWING  ADIT  ENTRANCE 
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FIGURE  9.  VIEW  OF  THE  COLORADO  SCHOOL  OF  MINES 

test  site  showing  the  rugged  surface  and  steep  slope 


FIGURE  10.  VIEW  OF  MINE  SHAFT  SHOWING  GRANITE  OUTCROPPING 


FIGURE  11.  CONTOUR  MAP  OF  THE  AREA  OVER 
THE  MINE  AT  IDAHO  SPRINGS,  COLORADO 


*6 


fLOOR  OF  MINE  IS  AT 
APPRO*  -<0 


(3  in)  of  soil  into  the  hole  leaving  a working  depth  of  50  feet  (15.  2 m). 

3.  Resistivity  Surveys  and  Results 

Three  resistivity  surveys  were  made  crossing 
the  mine  adit  where  overburden  thickness  was  30,  50  and  73  feet 
(9.  1,  15.  2 and  22.  2 m).  A set  of  measurements  was  also  made  down 
the  single  dry  borehole. 

The  pole-dipole  electrode  array  war  with 

both  AC-  and  DC-instrumentation.  Potential  electrode  spacings  of 
both  5 and  10  feet  (1.  5 and  3 m)  were  tried.  It  was  generally  found 
that  ten-foot  (3  m)  potential  electrode  spacing  gave  good  target 
resolution  but  the  five-foot  (1.  5 m)  spacing  caused  many  smaller 
anomalies  to  be  detected.  Most  t averses  were  run  using  a fixed 
potential  electrode  spacing  of  10  feet  (3  m)  incremented  m distance 
from  the  current  electrode  in  5-foot  (1.  5 m)  intervals.  The  results 
using  AC  and  DC  power  were  identical.  (This  comparison  was  made 
on  the  service  road  traverse  only.  ) 

Excellent  results  were  obtained  on  the  road 
traverse  where  the  adit  was  30  feet  (9.  1 m)  below  the  surface.  A total 
of  4 arc  pairs  crossed  in  the  volume  where  the  adit  was  located  as 
shown  in  Figure  13.  The  search  in  this  case  was  made  by  starting  with 
one  current  electrode  (infinity)  850  feet  (259.  1 m)  from  the  adit  center 
line  and  the  near  current  electrode  100  feet  (30.  5 m)  from  the  center 
line  in  the  opposite  direction  (refer  to  the  figure).  Starting  with  the 
10-foot-spaced  (3-m)  potential  electrodes  at  5 and  15  feet  (1.  5 and 
4.  6 m)  from  the  current  electrode  they  were  incremented  in  5-foot 
(1.  5-m)  intervals  to  the  10-  and  20-foot  (3-  and  6,  1-m)  marks  and  so 
on  out  to  the  65-,  75-foot  (19.  8-,  22.  9-m)  locations.  The  current 
electrode  was  then  moved  closer  to  the  adit  center  line  by  an 
increment  of  30  feet  (9.  1 m),  and  the  traverse  repeated  on  both  the 
forward  and  reverse  sides  of  the  current  electrode.  The  figure  shows 
the  locations  of  other  high  and  low  resistivity  manifestations  near  the 
road  surface.  The  low  resistivity  areas  labeled  "A"  and  "C"  were 
known  moist  areas.  The  high  resistivity  anomalies  were  caused  by 
near  surface  and  actual  outcroppings  of  granite. 

Construction  of  the  graphical  sketch  in  Figure  13 
was  done  as  was  shown  earlier  in  Figure  3 by  observing  the  graphs 
of  measured  data  for  high  and/or  low  resistivity  indications.  These 
were  then  labeled  and  their  locations  transferred  to  a graphical  sketch 
sheet  where  the  bearing  arcs  were  drawn. 
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Interpretation  of  the  data  taken  over  the  rough  broken 
rock  terrain  is  difficult.  Ground  surface  effects  are  very  strong,  and  in 
many  cases  they  almost  mask  the  anomaly  caused  by  the  mine  adit. 
However,  careful  observation  of  ground  surface  topography  aids  in 
eliminating  many  of  these  effects  from  the  data  graphs.  Figure  14 
shows  the  interpreted  results  of  a traverse  along  survey  line  A 
where  overburden  was  50  feet  (15.  2 m)  thick.  At  least  six  arcs  rep- 
resenting resistivity  highs  intersect  in  the  volume  of  rock  where  the 
mme  adit  is  located.  Also,  as  shown,  approximate  soil  depth  variations 
couid  be  determined  along  with  near  surface  rock.  Other  below  surface 
anomalies  were  not  verified,  but  the  low  resistivity  volume  seen  near 

che  mine  is  most  likely  caused  by  water.  There  was  seepage  into  the 
mine  in  this  area. 


The  results  of  a traverse  that  crossed  over  the 
mine  adit  where  it  was  73  feet  (22.  2 m)  below  the  surface  were  not  as 
outstanding.  Detection  results  of  the  traverse  along  survey  line  B 
are  shown  in  Figure  15.  There  were  high  resistivity  anomalies  noted 
at  locations  that  fitted  the  target  location;  but  a great  number  of  high 
resistivity  anomalies  (as  large  or  larger  than  the  mine  anomalies)  not 
associated  with  the  mine  adit  were  detected.  Most  of  these  "highs" 
were  surface  related  effects  from  observed  granite  outcropping.  On 
this  traverse  measurements  were  made  to  a distance  of  130  feet  (39  6 m 
each  side  of  the  current  electrode,  and  the  current  electrode  was  incre- 
mented in  40-foot  (12.  2 m)  intervals.  The  surface  anomalies  generally 
appear  on  the  graphs  of  three  consecutive  traverses,  have  about  the 
same  shape,  and  are  offset  from  curve  to  curve  by  a distance  of  10 
eet  (12.  2 m).  If  arcs  are  drawn  representing  these  "highs"  they 
will  intersect  at  the  ground  surface.  "Highs"  were  located  such 
that  there  were  arc  intersections  at  the  mine  location.  Near- surface 
resistivity  anomaly  interpretations  are  not  shown  in  the  Figure  15  sketch 

........  , . Resistivity  measurements  in  the  borehole  gave 

very  e mite  detection  indications.  Measurements  were  made  with  the 
current  source  electrode  first  at  the  top  of  the  hole  and  then  at  the 
bottom  of  the  hole.  Both  5-  and  10-foot  (1.  5-  and  3-m)  potential 
electrode  -pacings  were  used.  The  anomaly  appeared  larger  with 
the  5 -foot  (1.  5-m)  electrode  spacing. 

Data  interpretation  results  are  shown  in  Figure  16. 
ree  sectional  views  are  shown  so  that  any  effects  of  the  surface  slopes 
could  be  noted.  lnc  hole  was  drilled  normal  to  the  plane  of  the  road 
surface.  Because  of  the  slope  of  the  road,  the  borehole  had  an  angle 
that  was  13  from  vertical  and  also  angled  slightly  in  toward  the  mountain 
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BOREHOLE  DETECTION  RESULTS,  IDAHO  SPRINGS,  COLORADO 


IDAHO  SPRINGS,  COLORADO 


BOREHOLE 


FIGURE  16  (cont.  ) SECTIONAL  VIEWS  OF 
BOREHOLE  DETECTION  RESULTS,  IDAHO  SPR INGS,  COLORADO 


, i . Figure  16(a),  Section  A,  is  a plane  through  the 

locTtioIsr  ThI 1 rr  '1°  thC  SGrViCe  (SGe  Fi8Ure  12  f°r  Sectional 

surf_  ‘ 'S  Partlcular  section  provides  a view  of  the  maximum 
surface  slope  that  might  affect  resistivity  measurements.  In  this 
view  the  mine  adit  appears  horizontal  running  from  left  to  right  on  the 

uof„Ct'ed  tr“0nK  rrUUS  *”  Sh0Wn  in  thC  ««“«•  The  locations 
P in  ed  out  are  high  resistivity  arc  crossings  that  are  in  the  correct 

location  to  be  caused  by  the  mine.  The  two  marked  high  resistivity 

DOtCmtri<T1b8nat  thV1Rht  °'  the  Pa«e  are  sections  of  the  same 
potential  shells  as  those  shown  at  the  left  of  the  page.  It  must  be 

remembered  that  when  the  current  source  is  below  the  ground  surface 

s hZT  shells  around  it  are  spherical  in  shape  if  the  medium 

hrouThout  M°US'  Hom°80ncous  subsurface  conditions  were  assumed 
throughout  these  tests.  Detection  from  at  least  three  holes  is  neces- 
ar>  to  obtain  an  absolute  location  of  the  detected  target. 

a ,,  The  hiSh  resistivity  anomaly  caused  by  the  mine 

adit  was  partially  countered  by  the  low  resistivi*w  ^ nnmeal  • u 

Had  these  not  been  encountered  the  entire  volume  around  the  mine" 

tTvUv  ind^  ?PPCared  35  Mgh  resistivitV  indications.  The  low  resis- 
t>  indication  surrounding  the  borehole  at  the  30-  to  40-foot  (9  1- 

, - fepth  was  caused  W water.  Water  was  found  at  this  depth 

when  the  Im.e  was  drilled!  however,  i,  drained  before  ,he  meturemen.s 

j.  , , , Thls  sectional  view  shows  that  the  surface  slone 

s rCeVh:  hiehU  dC‘CCli°”-  Had  crossed  fhe 

The  , f 8 resistivity  indication  would  have  been  questionable 

borehn,  G °n  eqU1P°tential  arcs  centered  at  the  top  of  the 

apparent  ^ ~ng 

Ihe  borehole  and  parpe^u^  “ ^ 

::  :: t'™.*  t thc bor'ho10 to ,h* a<"*-  ■> »tr;r„:i; 

.ho  d'rcc“on  ,s  ‘e™  10  h»  relatively  flat  compared  to  that 

IrZdTh  Ure  " 'V:  ThC  r“U,W“V  "itth  anomalies  can  be  seen 

“ lhlS  VieW'  0the™i8-  “he  ‘»e  same  as 


•hrough  the  borehole  ‘ TgaiT 
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In  summary,  the  mine  adit  was  detected  from 
resistivity  measurements  made  in  a borehole  that  passed  the  adit  at 
a distance  of  16  feet  (4,  9 m).  Best  results  were  obtained  when  the 
potential  electrode  spacing  in  the  pole-dipole  electrode  array  was 
5 feet  ( 1 . 5 m). 


B, 


Gold  Hill,  Colorado,  Site 


1. 


Site  Description 


The  town  of  Gold  Hill,  Colorado,  is  located 
approximately  15  miles  (24  km)  northwest  of  Boulder,  The  Gold  Hill 
test  site  is  located  about  one  mile  (1.  6 km)  west  of  the  town  at  an 
elevation  of  8500  feet  (2.  6 km).  It  is  the  site  of  an  abandoned  mine  that 
has  an  adit  that  runs  through  the  mountain  a distance  of  approximately 
600  feet  (183  m).  The  cross  section  of  the  passage  is  about  5x6  feet 
(1.  5x1,8  m).  The  ground  surface  has  an  average  upward  slope  from 
the  adit  entrance  of  23°.  The  ground  surface  is  covered  moutly  with 
grass  and  large  pine  trees.  There  is  very  little  brush  to  interfere 
with  surface  surveys.  The  subsurface  structure  is  granite  overlain 
by  a sandy  mixture  of  mostly  weathered  granite.  A large  amount  of 
quartz  is  scattered  over  the  ground  surface.  Figure  17is  a view  of  the 
area  looking  southeast  along  survey  line  Number  Three  where  one 
resistivity  search  traverse  was  made.  The  mountain  in  the  background 
was  the  location  of  the  infinity  current  electrode. 

Two  faults  ran  through  the  test  area  almost  parallel 
and  about  150  feet  (45.  7 m)  apart.  Resistivity  survey  traverse  lines 
crossed  these  fault  lines  in  several  places.  Mining  had  been  done  along 
t ese  faults  leaving  craters  of  various  sizes  along  the  way.  They 
averaged  about  30  feet  (9.  1 m)  in  diameter  with  depths  o(  about  20  feet 
(6.  1 m)  to  the  partially  filled  in  bottom.  Figure  1 Sis  a view  of  a 
typical  crater. 

Comparing  the  Gold  Hill  test  site  to  the  Colorado 
School  of  Mines  site,  the  Gold  Hill  site  was  less  rugged,  and  slopes 
were  not  as  steep.  It  also  had  more  top  soil  making  electrode  place- 
ment easier.  The  many  craters  at  Gold  Hill  caused  resistivity  anomalies 
when  a traverse  came  near  one,  but  these  were  identifiable  when  the  data 

were  processed.  Overall,  the  Gold  Hill  site  was  more  accessible  and 
easier  to  work. 
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FIGURE  18.  VIEW  SHOWING  A LARGE  CRATER  NEAF 
RESISTIVITY  TRAVERSE  LINES  AT  THE  GOLD  HILL, 
COLORADO,  TEST  SITE 


2. 


t 


Site  Preparation 

The  only  site  preparation  necessary  was  to  drill 
holes  for  borehole  testing.  The  United  States  Bureau  of  Mines  in 
Denver,  Colorado,  had  an  interest  in  the  results  of  the  tests  to  be 
done  at  the  Gold  Hill  site.  They  had  used  the  site  for  some  of  their 
experiments  and  obtained  permission  for  work  on  this  program  to 
be  done  there.  Survey  lines  were  already  staked  from  their  previous 
work,  and  these  were  useable  for  the  resistivity  surveys.  The 
Bureau  of  Mines  al.o  needed  boreholes  for  their  tests.  They  performed 
the  drilling  and  al  owed  the  holes  to  be  used  for  tests  on  this  program. 

Two  boreholes  were  planned:  one  25  feet  (7.  6 m) 

from  the  adit  center  line  and  one  50  feet  (15.  2 m)  from  the  center  line. 
These  were  to  be  drilled  at  a location  where  there  was  about  80  feet 
(24.  4 m)  of  overburden  and  drilled  to  a depth  of  160  feet  (48.  8 m). 

After  a great  deal  of  equipment  trouble,  one  hole  3 inches  (76.  2 mm) 
in  diameter  was  completed  to  a depth  of  135  feet  (41. 2 m)  by  the  time  the 
field  measurements  had  to  be  completed.  Both  holes  were  finally 
completed  to  160  feet  (48.  8 m)  but  too  late  to  be  used  in  the  current 
program. 


3-  Resistivity  Surveys  and  Results 

Figure  19  is  a topographical  sketch  of  the  test 
are;  . The  adit  runs  from  the  road  level,  Hazel  A Portal,  into  the 
mountain  (bottom  of  the  page  to  the  top).  The  Bureau  of  Mines  survey 
lines  numbered  3,  4,  and  7 are  shown  crossing  the  adit.  The  borehole 
locations  are  shown  on  Line  7.  Overburden  under  Line  3 was  30  feet 
(9.  1 m);  under  Line  4,  47  feet  (14.  3 m);  and  under  Line  7,  81  feet 
(24.  7 m).  Resistivity  traverses  were  made  along  these  lines  and 
down  the  one  partially  completed  borehole. 


Using  potential  electrode  spacing  of  5 feet  (1.  5 m) 
incremented  in  5-foot  (1.  5-m)  intervals  the  mine  adit  was  detected  on  all 
three  surface  traverses.  Not  only  was  the  adit  detected,  but  all  known 
surface  features  could  be  distinguished  by  their  resistivity  anomalies. 

An  example  of  detection  results  using  the  graphical  construction  method 
is  shown  in  Figure  20.  The  traverse  shown  is  along  Line  7 on  the 
topographic  map  (see  Figure  19  ).  The  map  shows  two  rather  large 
patches  of  granite  outcropping  starting  about  60  feet  (18.  3 m)  to  the 
right  of  the  adit  center  line.  Also,  two  fault  lines  cross  the  traverse. 

As  can  be  seen  on  the  resistivity  data  sketch,  Figure  20,  there  are  large 
concentrations  of  high  resistivity  indications  at  both  granite  outcrop 
locations.  Also  noted  are  strong  low  resistivity  indications  where  the 
faults  were  crossed  (arcs  not  drawn  in  purposely  to  reduce  clutter). 
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SCALE  : 1 in  - 50ft. 


3 


FIGURE  19.  TOPOGRAPHICAL  MAP  OF  THE 
GOLD  HILL,  COLORADO,  TEST  AREA 
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The  "lows"  were  probably  caused  by  moisture  concentration  and  other 
conducting  materials  in  the  faults. 

The  mine  adit  showed  up  very  clearly  even  at 
the  81-foot  (24.  7-m)  depth.  Six  arc  pairs  crossed  at  the  adit  location. 

An  interesting  detection  result  is  the  detection  of  the  subsurface  water 
in  the  volume  of  earth  where  the  hole  uas  bored.  This  was  a very 
positive  indication  with  8 arc  pairs  crossing  in  that  location. 

Figured  21  and  22  are  included  to  provide  an  example 
of  the  data  plots  from  which  the  interpretation  in  Figure  20  was  made.  On 
each  sheet  the  current  source  electrode  positions  are  marked  on  the 
horizontal  line  labeled  SURFACE.  The  curves  immediately  to  the  left 
and  right  of  each  current  electrode  position  are  the  data  curves  plotted 
from  traverses  run  to  the  left  and  right  of  that  current  electrode 
position.  The  traverses  on  which  these  data  were  taken  were  nominally 
200  feet  (6  l m)  on  each  side  of  the  current  electrode.  In  some  cases 
when  the  traverse  was  directed  away  from  the  adit,  the  traverse  was 
shorter.  In  other  cases  the  traverse  length  was  extended  to  ascertain 
that  it  crossed  the  adit  far  enough  to  allow  detection. 

Only  the  anomalies  caused  by  the  adit  are  labeled. 
These  were  determined  when  the  data  were  studied  and  analyzed.  The 
small  size  of  the  mine  anomaly  compared  to  other  anomalies  will  be 
noted.  On  analyzing  the  data,  it  will  be  found  that  most  of  the  large 
anomalies  are  caused  by  near  surface  rock  outcropping  or  by  craters 
near  the  traverse  line.  It  must  be  remembered  that  the  potential 
bowls  or  resistivity  shells  around  the  current  electrode  are  shells 
somewhat  hemispherical  in  shape  (hemispherical  in  a homogeneous 
medium),  and  resistivity  anomalies  within  several  feet  of  a traverse 
line  will  be  detected. 


Detection  of  the  adit  under  30  feet  (9.  1 m)  of 
overburden  was  definite.  At  the  47-foot  (14.  3-m)  overburden 
location,  detection  indications  were  strong;  but  a crater  only  a few  feet 
southwest  of  the  adit  center  line  (refer  to  Figure  ’9)  on  this  traverse 
confused  the  results,  because  data  arcs  p?  ssing  through  the  adit  also 
passed  through  the  crater.  Anyone  interpreting  the  data  from  that 
location  not  knowing  the  adit  was  below  would  have  misinterpreted  the 
results  as  being  caused  by  the  crater. 

Borehole  resistivity  measurements  were  mostly 
unsuccessful.  The  borehole  had  hit  water  at  a depth  of  40  feet  (12.  2 m) 
and  remained  full  to  that  level.  It  was  pumped  dry  but  continued  to 
refill  too  fast  to  allow  measurements  to  be  made  in  a dry  hole.  Since 
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FIGURE  21.  RESISTIVITY  DATA  GRAPHS  OF  THE  SOUTHWEST  END 
OF  LINE  7 TRAVERSE,  GOLD  HILL,  COLORADO 


SURFACE 


FIGURE  22.  RESISTIVITY  DATA  GRAPHS  OF  THE  NORTH  EAST  END 
OF  LINE  7 TRAVERSE.  '''•OLD  HILL,  COLORADO 


““  be  kept  dry’  “ Was  kept  fu"  o!  water  to  the  ground  surface 
so  that  all  measurements  would  be  made  in  water  Th!  7 

I"  drilling  this  hole  th^Ml',^'1’  *ekWalU  °f  tha  hale  was  also  a problem, 
all  drill  rods  tnat  w erl  put  in  o th  't  , r"  °f  thick  CUP  Urease  on 

drill  rod  and  make  drilling  easier  Tts  ? ^ ”U  MM°"  ‘ho 
that  has  a high  electrical  re^r  * Th  ® left  a very  heavy  coat  of  grease 

one  indication  of  detection  ^ ^ b°reh°'e  showed 

other  overlapping dlta  o actuall  WaS  P°SltWa  b'“"se  there  was  no 
aaca  to  actually  prove  detecHnn  ti 

resistivity  high  that  occurred  at  the  cor rect  1 l“„ce  from  7b  “ °n<! 
electrode.  aistance  from  the  current 


additional  trip  to  the  si^inTThi  Tl  aVaUaMe  did  Permit  an 

surveys  in  water-filled  holTc^d ^ZVn 
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IV-  CONCLUSIONS 

experimentally^  eva hia t’e  larth^e^T-  dlSCUSSed.in  this  report  was  to 
subsurface  voids  in  a granite  environmlm^  Thl"15  * melhod  of  detecting 
a high-resolution  electrode  arrav  It-  ' resistlvity  method  used 

electrode  array  AUhlugh  " “ 8 P°le-diP°la  »'  three- 

cone, usions  were  reached as  ^ "*  8 "a™b«  »' 

d-8  x 1.8  m)  can  be  detected^'t  do7h  ” .SeCti°ns  in  the  range  6x6  feet 
over  80  feet  (24  4 nil  The  m • P ^ be  tho  ground  surface  of 
probably  about  100  feet  ,30  5 m^T,  n P‘h  WaS  "°l  dotermined.  La, 
tunnels  with  the  e,e‘v^“ LT  ^ 

does  no,  appea^to^ffecr det^^tio^csipabiUtie^s1.^  S'°P°S  °‘  ^ 3°‘ 
variable  soil  d/pXTulf.TaVd^^^^ 

pretatxon  very  complex-  hut  in  6 , can  make  data  xnter- 

method  allows  these  features  to  be  mapped!  C reS°1Uti°n  °f  the  Section 


4. 

variations. 


The  method  is  most  sensitive  to  near-surface  resistivity 


sensitive  the  mltt^'e^m^VC^'  Sp8ai"S  the  less 

smaller  the  spacing  with  respect  to  th^H*-  varlatlons;  aIso*  the 

the  larger  the  detection  indication.  TherefTre  th°f  Subsurface  larget 
electrode  spacing  would  be  la  rao  , 1 e’  the  optimum  potential 

variation  (noise)  as  possible  but  smaUeno  ^hT^u  ^ background 

capability.  A spacinPg  no  sLlllrThan  7 ° Cnhance  Urget  detectio- 

suspected  target  should  be  an  approprUteecc^pr^i7edian,eteP  * ^ 

array  as  described  in  this  repo^^is  ^ three_electrode  resistivity 
two-man  crew  can  complete  a single  ZOQ^n  f°™UminZ  Pr°cess.  A 
about  IS  minutes;  however  the  rem  i h j -m)  survey  traverse  in 
measurement  will  increase’  the  ^ ^ °‘ 

average  one  500-foot  (152.  4 m)  survey  traverse’  in  an“r  daT 


Preceding  page  blank 
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.,  . Drilling  is  the  only  known  method  to  'tbsolutely  determine 

the  existence  of  air-  or  solution-filled  subsurface  cavities.  Resistivity 
measurements  can  greatly  minimize  the  required  numbex  of  boreholes 
in  an  area  by  identifying  anomalous  areas  that  could  be  cavities  or  tunnels. 


The  three  electrode  resistivity  measurement  method  of 
cavity  detection  can  be  used  in  boreholes.  Although  this  arrangement 
ad  a very  limited  test  on  this  program,  a mine  adit  having  a cross 
sectmn  of  about  10  x 10  feet  (3  x 3 m)  was  detected  from  a dry  borehole 
16  feet  4.  9 m)  away.  It  is  probable  that  such  detection  could  be  made 
frorri  hole  distances  of  up  to  40  feet  (12.  2 m).  This  use  shows  promise 
but  for  better  understanding  and  detection  range  prediction,  both  theoret- 
ical and  experimental  studies  should  be  continued. 


, rpun,The,n  lre  tW°,  great  advanta8es  to  making  measurements  in  a 
borehole:  (1)  the  rock  environment  near  the  electrodes  should  be  much 

fOUnd  near  the  gr°Und  SUrface*  thus  eliminating 

indi  r m f 81Cal  n°ise"  and  enhancin8  the  desired  anomaly 

ndication,  (2)  if  the  target  tunnel  is  very  deep,  a borehole  can  put  the 

e ectrodes  within  the  detection  range  of  the  search  method. 


V. 


recommendations 


cavitiesTa„dtthurnne',eltCtr°Kde  rC8isti'’il''  method  °'  detecting  underground 

It  has  detef  *ad  H Va  proven  successful  in  several  environments. 

through  sandv  lV°  , P 5 °f  over  80  '<==■  * ™>  and  has  worked 

rough  sandy  loam  soil  over  limestone,  sand  over  limestone,  moist 

y over  dolomite,  and  through  granite.  Sometimes  detection  was 

to'Imoro  d"f  °!her  timCS  “ was  marginal.  Several  recommendations 

can  be  made  b-ed  "■•••  «•« 


A. 


Additional  Experiments  and  Studies 


for  rh»  I it  • Speclflc  experiments  and  analyses  should  be  performed 

for  the  following  purposes:  (a)  to  determine  the  ultimate  spatial 

resolution  the,  can  be  obtained;  ,b)  to  determine  the  effect  v volume 

a SUV'  ' rd  fieM  °f  WeW  affeCti"8  measurements^  .1™* 

a survey  path;  (c)  to  improve  or  develop  better  methods  of  field  data 

n erpretation;  and  (d)  to  study  other  electrode  arrangements  that 
might  enhance  tunnel  anomaly  indications.  8 that 

It  was  stated  by  both  Bates4  and  Bristow5  that  when 
usmg  he  pole  -dipole  electrode  array  in  the  manner  reported 'herein 
e volume  of  earth  material  affecting  the  potential  measurements  ’ 

tt  a„n.t“C  O I"  ^ tH:  materUl  impHed  ^ thG  “ ed  in 

sheiu  , y a SeCt°r  0f  the  graphical  hemispherical  potential 

I?ne  T *Z  than  ab°Ut  ± “ d^rees  each  side  of  the  traverse 

tMs  progral  h l°  the  Measurement.  Results  obtained  on 

his  program  have  roughly  validated  this  limit  although  no  detailed 

interpreTaaLTV1316'  ^ 3 Very  imPortant  operational  and 

of  resisUvitv  me  * " Sh°Uld  be  Studied  bV  laboratory  modeling 
b.  additional  fieM?^”16"15  Z 3 brine  Unk  °r  °ther  suitable  medium, 
analysis  The  sa  ? °Ver  SubsUrface  targets,  and  by  theoretica 

spatial  resolution  a^^^^ 


Rates,  loc.  cit. 


Bristow,  loc.  cit. 
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to  improved  graphical  analysis  and,  possibly,  by  computer  processing, 
to  automatic  data  reduction  and  subsurface  mapping. 

Other  field  data  collection  methods  might  be  used  that 
would  take  advantage  of  the  physical  shape  of  a tunnel.  Since  the  tunnel 
can  be  visualized  as  a long  cylinder,  an  electrode  array  could  probably 
be  used  that  would  take  advantage  of  tnat  shape.  Line  arrays  could  be 
laid  out  and  measurements  made  in  the  same  manner  as  was  described 
for  the  standard  pole-dipole  meas urement.  For  example,  a current 
electrode  in  the  form  of  a wire  100  feet  (30.  5 m)  long  could  be  laid  out 
parallel  to  a suspected  tunnel  axis  and  staked  to  the  ground  at  several 
points.  Potential  electrodes  of  the  same  form  would  then  be  used  to 
make  potential  measurements.  Such  an  array  should  have  a signal  averaging 
effect  that  would  enhance  the  tunnel  anomaly  and  reduce  near  surface 
clutter  not  common  to  the  entire  length  of  the  line  electrodes. 

Other  electrode  arrays  such  as  a dipole -dipole  array 
should  also  be  investigated. 

B.  Study  the  Resistivity  Contrasts  Between  Earth 

Materials  Associated  With  Cavity  Detection 

The  size  of  the  resistivity  anomaly  that  a given  subsurface 
target  will  cause  can  only  be  derived  if  the  resistivities  of  the  various 
materials  are  known.  Samples  of  soil,  bedrock,  and  water-soil  mixtures 
from  water-  and  mud-filled  cavities  should  be  collected  from  areas  of 
interest  and  their  electrical  and  electromagnetic  properties  measured 
and  analyzed.  A great  deal  of  information  is  available  on  electrical 
properties  of  rock,  but  more  data  are  needed  on  those  materials 
covering  and  in  the  vicinity  of  the  rock. 

C-  Develop  State-of-the-Art  AC  Resistivity  Inst rumentation 

Ground  current  variations  during  the  course  of  a survey 
traverse  can  cause  errors  that  appear  to  be  resistivity  anomalies  when 
the  data  are  analyzed.  A temperature  stable  constant-current  power 
source  should  be  developed  to  work  in  a wide  range  of  ambient 
temperatures  and  that  can  supply  ground  currents  under  very  dry  soil 
conditions.  It  should  be  capable  of  supplying  a constant  current  ot  at 
least  10  ma  when  combined  earth  resistance  and  electrode  contact 
resistance  is  as  high  as  50,000  ohms. 
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. v°ltnioter  section  of  the  instrumentation  should 

>e  built  to  have  a very  high  input  impedance,  in  the  range  of  JO  megohms, 
to  prevent  potential  measurement  errors  that  might  be  caused  by  high 
electrode  contact  resistance.  A low  impedance  meter  across  the  potential 
electrodes  would  allow  current  to  flow  in  a circuit  through  the  meter 
the  electrode  contact  resistance,  and  the  earth,  causing  a potential  drop 
across  any  electrode  contact  resistance.  This  would  produce  an  error 
voltage  which  might  be  interpreted  as  measured  earth  potential.  Such 

eatures  as  narrow-band  synchronous  detection  should  be  built  into  the 
meter  input  circuit  to  give  best  possible  signal-to -nois e ratio.  To  aid 
in  making  field  measurements  as  fast  as  possible,  an  autoranging  volt- 
meter  should  be  used  with  a digital  readout.  Another  feature  that  should 
be  added  especially  for  borehole  operation  is  a recorder  output  so  that 
continuous  data  could  be  recorded  as  electrodes  were  moved  up  and  down 
the  borehole  walls. 


D. 


Develop  a Mobile.  Rolling- Contact- Electrode 
Earth  Resistivity  Survey  System 


. , 11  1S  recornniended  that  an  earth  resistivity  system 

be  developed  that  would  allow  continuous  measurements  using  a pole- 
dipole  electrode  arrangement  in  which  the  electrodes  are  pulled  along 
the  ground  surface  using  a small  tractor-type  vehicle.  Automatic  data 
acquisition  and  processing  could  readily  be  built  into  such  a system  to 

variatior^s  re<H  ^ra^11Cal  printou^s  or  other  displays  of  apparent  resistivity 
variations  and  depth  profiles  as  the  electrodes  move  along  a selected 

/averse.  With  digital  data  processing  the  graphical  analysis  as  described 

and  TT  reduction  could  provide  near  real  time  output  information 

hardcopy  maps  delineating  the  locations,  depths  and  sizes  of  the 
detected  underground  cavities. 


Another  version  of  continuous  measurements  for 
errain  where  a vehicle  cannot  travel  would  be  a system  whereby  current 

tvnes  nS-arC  P , ^ USUal’  ^ thC  P°tential  electrodes  are  rolling 

th at  a rom™  ^ ! ***  °f  ab°Ut  5 ^ 5 m)  lt  is  conceivable 

f O 1H  h > array  comPrisin8  °"ly  the  two  potential  electrodes 

a ma  \ 1 ‘ ^ PUShed  °r  puUed  alonR  the  ground  surface  by 

a man  where  he  terrain  would  allow.  The  output  from  the  array  could 

be  continuously  recorded  for  later  interpretation.  The  ability  to  measure 

r h resistivity  using  rolling  electrodes  was  demonstrated  in  1 972  by  SwR] 


L.  S.  Fountain,  "An  Exploratory  Study  of  Soil  Resist 

sing  a Rolling  Contact  Electrode  Array,"  Final  Tech 
Internal  Research  Project  14-9057,  April  1972. 


ivity  Measurements 
nical  Report,  SwRI 
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E.  Continue  the  Investigation  of  Resistivity 
Borehole  Techniques 

1 his  would  require  the  development  of  better  borehole 
electrodes,  collection  of  much  more  field  data  in  areas  of  known  targets 
such  as  mine  adits  or  caves,  and  solving  the  problem  of  making  measure 
ments  in  a liquid-filled  hole. 


E.  Lndertake  the  Long-Range  Development  of  a State-of- 
the-Art  Ground  Penetrating  Radar  for  Tunnel  Search 
Applications 

, Radar  tests  have  shown  very  good  possibilities  for  the  method  7 
The  main  problems  encountered  in  the  field  have  been  penetration  depth 
and  small  target  detection  and  resolution.  Radar  technology  is  presently 
available  to  permit  the  development  of  a ground  penetrating  radar  system 
m ving  Significantly  improved  performance  over  any  that  are  presently 
in  use.  The  detection  depth  capabilities  could  be  expanded  greatly  by 
increasing  transmitter  power,  by  improving  receiver  design,  and  by 
incorporating  advanced  radar  signal  processing  techniques. 


While  the  pole-dipole  earth  resistivity  method  has 
demonstrated  an  early  and  successful  capability  for  detecting  subsurface 
cavities  and  tunnels,  the  ground  penetrating  radar  technique  holds  a 
ar  greater  potential  for  providing  more  detailed  subsurface  resolution 
and  a faster  ground  scanning  rate.  The  development  of  a ground  penetrating 
'adar  system  capable  of  reliably  scanning  and  accurately  mapping  under- 
ground cavities  and  other  subsurface  anomalies  is  recommended  as  a 
most  favorable  long  range  approach  to  solving  problems  related  to 
subsurface  cavity  detection. 


, A rnodular  radar  system  design  is  recommended,  allowing 

sUectlon  of  slPnal  frequency  and  matching  antenna  to  best  match 
the  system  to  the  soil  and  rock  medium  and  penetration  depth  required 
in  the  selected  urea.  This  would  include  special  antenna  development  for 

bo  rehole  ***”  ^ PCnetrate  thrc^h  thc  walls  of  a 

ehote  for  a distance  range  of  30  feet  (9.  J m)  or  more  would  be  of  great 

a ue  when  suspected  tunnels  are  too  far  below  ground  surface  to  be 
detected  by  surface  methods. 


Fountain.  Herzig,  Owen.  loc.  cit. 
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BASIC  CONCEPTS  OF  ELECTRICAL  EARTH 
RESISTIVITY  MEASUREMENTS 


The  basic  resistivity  measurements  and  analysis  concepts  for 
the  venous  survey  methods  are  essentially  the  same.  If  an  otherwise 
predictable  electrical  current  distribution  is  established  within  a rela- 
tively large  volume  of  homogeneous  earth  material,  perturbations  in  the 
current  distributions  as  caused  by  localized  geologic  anomalies  can  be 
measured  as  potentials  at  the  ground  surface.  These  anomalous  poten- 
tials can  then  be  interpreted  in  terms  of  possible  subsurface  structural 
or  earth  materials  differences. 

The  basic  electrical  theory  for  earth  resistivity  measurements 
is  discussed  and  the  geometrical  factors  for  several  electrode  arrays 
are  derived  below.  y 


I. 


GENERAL  THEORY 


An  electric  current  distribution  can  be  established  within  the 
volume  of  a conducting  earth  haxfspace  by  means  of  a DC  or  low  fre- 
quency  AC  voltage  applied  to  two  electrodes  in  contact  with  the  surface 
If  the  two  electrodes  are  located  close  together  the  volume  of  earth 
material  carrying  the  major  portion  of  the  total  current  will  be  relative- 

^diamete^ nd.Lanf  bC  ^mated  as  a hemispherical  volume  having 

meter  about  equal  to  the  electrode  spacing  if  the  composite  earth 
a enal  can  be  considered  conductively  homogeneous.  For  larger 
! °de  separation  distances  the  depths  of  the  significant  current  con- 
ducting paths  increase  as  does  the  volume  of  material  carrying  the 

^°r^T10n,0f!>  e t0taluCUrrent  Under  the  idealized  condition  of  a 
umforrtUy  conducting  earth  halfspace  and  in  the  limit  where  the  separa- 
tion distance  between  the  two  electrodes  approaches  infinity  (i  e spac 
lng  '' e ry  lar«e  compared  with  the  dimensions  of  the  earth  material  zone 
containing  the  current  distribution  patterns  of  interest),  the  current 
owing  away  from  a point  source  electrode  diffuses  radially  into  the 

liderelT  ^ ideali‘ed  case’  the  current  flow  can  be  con- 

sidered as  diverging  away  from  the  source  through  an  infinite  number  of 

“ ^ hTSPheTal  SheUS  “P™"***  equipotential  surfaces 
n the  medium.  At  the  ground  surface,  the  hemispherical  sheUs 
orm  concentric  circles  about  the  current  electrode  at  which  potential 
Inferences  between  such  equipotential  shells  may  be  measured. 

By  Ohm's  law,  the  vector  potential  gradient  between  two  infini- 
tesimally spaced  hemispherical  shells  is 


Preceding  page  blank 
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(A-l) 


where  p = resistivity  of  the  medium;  and 

j = vector  current  density. 


For  a single  current  source  with  its  current  sink  at  infinity,  all  current 
entering  the  hemispherical  volume  from  the  source  must  also  emerge 
from  that  volume  and,  therefore,  the  vector  current  density  divergence 
condition 


V • 


(A -2) 


must  be  satisfied.  Substituting  from  Equation  (A-l),  the  Laplace  equa- 
tion for  the  scalar  electric  potential  is  found  to  be 

7 ' T - j v - f,0 

or 

(A -3) 

— v-2  0 = 0 

r — J3 

where.  <*>  =J  £ • d f = scalar  electric  potential. 

OO 

On  the  basis  of  the  radial  (hemispherical)  symmetry  already  described 
for  a homogeneous  conducting  earth,  Laplace's  equation  in  spherical 
coordinates  is  simply 


d_ 

dr 


(A -4) 


which,  after  direct  integration,  yields  the  electric  field  relation 


r 


£ = x, 


(A-5) 


and, 


after  integrating  again,  the  electric  potential  relation 


0 = - 


+ k2 


(A  -6) 
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where:  Kl,  K2  = constants  of  integration. 

At  large  distances  away  from  the  current  source  electrode  (r-*-  °°  ) the 
potential  0 must  be  zero  and,  therefore,  K2 — ►O  . Further,  if  the 
source  current  is  uniform  through  a small  hemispherical  surface  about 
the  source  /’ectrode  then  the  total  current  is  given  by  the  surface 
integral 


1 "/  l f C * ds  • 

C ' n 


(A-7) 


Substituting  from  Equation  (A -5),  the  total  current  becomes 


Ki  f I Ki/*  f 2ttP 

I = f — ds  = J dy  J sin  6 dO  

P { r*  P o -rr/2  P 


(A -8) 


where  ds  = r^  sin  0 d-ydO; 

y - azimuthal  angle  about  current  source;  and 
0 = polar  angle  about  current  source. 

Solving  for  Kx  and  substituting  into  Equation  (A-6)  gives  the  scalar  elec- 
tric potential  for  any  hemispherical  shell  of  radius,  r,  about  a single 
current  source  carrying  a total  current,  I.  That  is, 


0 f r ) = . Rl 

2 tt  r 


(A -9) 


The  observable  potential  difference  between  two  points,  Pj  and  P2,  on 
hemispherical  shells  concentric  about  the  current  source  and  located  at 
radii,  r j and  r2  where  r2  > r j , is 


A 0 = 0 (rj ) - 0 (r?)  = -M—  - - _L 

L rj  r 


(A -10) 


This  relationship  is  important  in  defining  the  resistivity  of  the  earth 
medium  *hen  the  current  sink  electrode  is  located  effectively  at  an 

infinite  distance  away  from  the  current  source  electrode.  In  this  case 
the  resistivity  is  ’ 


= K (■ 


(A  -1 1 ) 
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= electrode  array  geometric  factor;  and 


^0  = a measurable  ohmic  resistance  factor. 

I 

For  the  case  where  the  current  sink  electrode  cannot  be  consid- 
ered to  be  located  at  an  infinite  distance  away  from  the  source  electrode, 
the  two  equipotential  shells  used  in  obtaining  Equation  (A -10)  will  no 
longer  be  precise  hemispheres  and  the  distortion  effects  on  the  poten- 
tials 0 ( r i ) and  0(r2>  caused  by  the  sink  electrode  must  be  taken  into 
account.  This  can  be  done  readily  because  of  the  validity  of  superim- 
posing scalar  quantities  such  as  the  electric  potentials  as  expressed  by 
Equation  (A -9).  Thus,  if  the  total  sink  electrode  current  is  -I  as  it 
must  be  in  a two -elect rode  arrangement  and  the  distances  from  the  sink 
electrode  to  the  previously  described  potentialmeasuring  points  (located 
at  distances  of  rj  and  r2,  respectively,  from  the  source)areRj  andR;?, 
respectively,  then  the  combined  potentials  at  the  two  measuring  points 
are  (for  Rj,  R2  > rj,  r^), 

0 pj  = 0 f r i ) - 0 fR  i ) 

and  (A -12) 

0 pz  = 0(r2)  - 0 (R2)  ■ 

The  resulting  potential  difference  between  the  two  potential  measuring 
points  now  becomes, 


£0  = 


r 2 + R 2 


(A -1 3) 


and  in  this  more  generalized  case  the  resistivity  of  the  medium  is  ex 
pressed  as 


P 


2u 

J_  - _L  - J_  + _L 

r j R,  R2 


).k 


(A -14) 


where  the  geometric  factor,  K',  for  this  case  differs  from  K expressed 
earlier  in  Equation  ( A - 1 1 ) to  take  into  account  the  current  sink  elec- 
trode distances. 


It  is  pointed  out  that  the  relationships  given  in  Equations  (A -13) 
and  (A-14)  are  general  expressions  for  any  electrical  resistivity  pro- 
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filing  array  which  utilizes  two  current  electrodes  and  two  potential 
electrodes  as  long  as  the  distances  rj  and  r2  represent  the  magnitudes 
of  the  geometric  distances  from  potential  electrodes,  P],  P2,  respec- 
tively, to  the  current  source  electrode  and  distances  Rj  and  R2 
represent  the  magnitudes  of  the  distances  from  potential  electrodes, 
Pl,P2»  respectively,  to  the  current  sink  electrode,  C2- 


II.  THE  WENNER  ELECTRODE  ARRAY 

The  Wenner  electrode  array  is  one  of  the  most  widely  used 
arrays  for  measuring  earth  resistivity.  This  array  is  a four-electrode 
configuration  in  which  all  electrodes  are  equally  spaced  along  a straight 
line.  The  distance  between  any  two  adjacent  electrodes  is  defined  as  the 
array  spacing.  The  Wenner  array  electrode  configuration  is  illustrated 
in  Figure  A-l.  The  geometrical  factor  for  the  Wenner  array  is,  from 
Equation  (A-14), 


Kw= = 2„d  (A-15> 

J_  _ J_  J_  + i_ 

d 2d  " 2d  "d 
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FIGURE  A-l.  THE  WENNER  ELECTRODE  ARRAY 
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The  geometrical  symmetry  of  this  array  makes  it  simple  to  set 
up  and  opeiate  in  the  field  and  the  simple  geometric  factor  facilitates 
the  interpretation  of  the  observed  electrical  resistivity  measurements. 

Referring  to  Figure  A-l,  if  the  current  sink  electrode  were  re- 
moved to  infinity,  the  hemispherical  equipotential  shells  associated 
with  potential  electrodes  Pj  and  P2  would  have  subsurface  penetration 
radii  of  d and  2d,  respectively.  That  is,  the  subsurface  earth  material 
affecting  the  measured  resistivity  is  that  contained  within  the  volume 
bounded  by  the  two  hemispherical  shells  concentrically  centered  on  the 
source  electrode,  C\.  However,  the  close  and  symmetrical  proximity 
of  the  current  sink  electrode  in  the  standard  Wenner  array  distorts  the 
equipotential  surfaces  away  from  their  concentric  hemispherical  shapes 
and  diminishes  the  depths  to  which  these  surfaces  penetrate.  This 
apparent  reduction  in  penetration  depth  for  the  standard  Wenner  array 
compared  with  the  single  current  electrode  case,  while  phenomenolo- 
gically correct,  is,  in  part,  compensated  by  the  larger  potential  dif- 
ference between  electrodes  Pj  and  P2  resulting  from  the  contribution 
of  the  current  sink  electrode.  The  apparent  maximum  penetration  re- 
sponse of  the  Wenner  array  is  oriented  directly  below  the  center  of  the 
symmetrical  array  configuration  in  contrast  with  the  location  below  the 
source  current  electrode  when  the  sink  electrode  is  located  at  infinity. 
In  practice,  the  standard  Wenner  array  has  its  optimum  response  to 
large  resistivity  anomalies  when  they  are  located  at  a subsurface  depth 
approximately  equal  to  the  array  spacing  dimf  nsion.  Because  of  the 
relatively  large  volume  of  earth  material  embraced  by  the  equipotential 
surfaces  terminating  at  electrodes  Pj  and  P2  and  surrounding  the  tv  o 
current  electrodes,  the  Wenner  array  tends  to  provide  a larger  volume 
resistivity  average  and  a lower  resolving  power  for  small  anomalier 
than  certain  other  electrode  configurations.  For  this  reason,  the 
Wenner  array  has  been  mi  st  applicable  in  geophysical  surveys  requiring 
depth  sounding  and  resistivity  measurements  related  to  large  subsurface 
earth  stratification  conditions.  While  it  has  also  been  used  as  a fixed- 
depth  horizontal  profiling  technique  (i.  e.  , array  spacing  held  constant 
and  the  e itire  array  moved  horizontally  over  the  ground  surface)  for 
detecting  vertical  earth  structural  formations  such  as  dikes  and  cavities , 
there  are  other  electrode  array  configurations  that  exhibit  g r eater  depth 
sensitivity  and  resolution  for  such  anomalies. 

III.  THE  POLE-DIPOLE  ARRAY 

The  pole-dipole  earth  resistivity  electrode  configuration  consists 
of  a current  source  electrode  (pole)  and  a potential  electrode  pair  (di- 
pole) oriented  in  a straight  line  array.  The  current  sink  electrode  is 
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located  sufficiently  far  from  the  source  electrode  that  its  effects  on  the 
potential  values  observed  by  the  dipole  electrode  pair  are  negligible. 

The  details  of  this  electrode  array  configuration  and  a description  of 
its  operation  is  presented  in  Section  IV  and  illustrated  in  Figure  9 in  the 
main  body  of  this  report. 

The  geometrical  factor  for  the  pole-dipole  array  is  derived  from 
the  general  expression  in  Equation  (A-14)  for  the  conditions  R1  — oo 
and  R2 -»  00  corresponding  to  the  location  of  the  current  sink  electrode 
at  an  infinite  distance  from  the  potential  electrodes.  Thus,  the  geomet- 
ric factor  for  the  pole-dipole  array  is, 


Kp-D  = 


Ztt 


1 . 


r2 


Zttti  rz 

r2  " r! 


(A -16) 


IV. 


THE  L-SHAPED  ARRAY 


The  L-shaped  electrode  array  is  comprised  of  a wide-spaced 
current  electrode  pair  located  along  the  lower  arm  of  the  "L"  pattern 
an  a closely  spaced  fixed -dimensional  potential  electrode  pair  oriented 
and  scanned  along  the  upper  perpendicular  arm  of  the  "L"  at  logarith- 
mically-spaced intervals.  The  source  current  electrode  is  located 
closest  to  the  potential  electrode  pair  at  the  corner  of  the  "L"  and  the 

the  •'L''CtTheiiS  ‘°Cat!d  a‘  the  °Uter  eXtr'me  point  °”  the  1°»"  a™  of 
L . The  L-shaped  an  ay  is  illustrated  in  plan  view  in  Figure  A -2. 

j , **  eeneral  the  current  sink  electrode  of  this  array  is  not  ronsid- 
ed  to  be  located  at  an  infinite  distance  away  from  the  current  source 
or  potential  measurement  electrodes.  However,  if  such  a condition  is 
pe  nutted  then  the  electrode  configuration  of  the  L-shaped  array  reduces 
other  ^ ^ already  described  for  the  pole-dipole  array  (and  for 

used  h array  COnfi^Urations  ^ the  half-Schlumberger  array 

used  by  previous  investigators).  The  exact  geometric  factor  for  the 
L-shaped  array  is,  from  Equation  (A-14) 


KL  = 


rl 


*1  r2 


/ 


r2  + b‘ 


2 +b2 


L j*!6  + /' 


(A  -1 7) 


r22  + b^ 
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1 


the  l-shaped  elelectrode  array 

FIGURE  A-2.  THE  L-SHAPED  ELECTRODE  ARRAY 


where:  b = distance  between  current  electrodes. 

As  may  be  noted  from  Equation  (A -17)  above  when  b — °°  the  geometric 
factor  for  the  L-shaped  array  approaches  that  given  earlier  in  Equation 
(A- 16)  for  the  pole-dipole  array. 

It  has  been  shown  by  Zohdy*  8 ) that  the  L-shaped  array  configu- 
ration has  advantages  and  flexibility  in  the  positioning  of  the  current 
sink  electrode,  especially  on  long  or  deep  sounding  resistivity  surveys, 
and  that  its  resistivity  sensitivity  to  stratified  earth  structures  is  inter- 
mediate between  that  of  an  ideal  equatorial  dipole  array  and  an  ideal 
polar  dipole  array  . That  is,  the  perpendicular  orientation  of  the  cur- 

A.  A.R  Zohdy,  "Electrical  Resistivity  Sounding  with  a L-Shaped 
Array,  U.  S.  Geological  Survey  B ulletin  1 3 1 3 - C,  1970. 


b °-le0  e?uatorial  array  is  one  in  which  the  distance 

dal  elecfr  dlStanCC  (r*  * r!>  ^ 0 *»  figure  A-2  and  then  the  poten- 

tial electrodes  are  repositioned  to  be  parallel  to  b.  , i.  e.  at  P,  anH 

P2-  An  ideal  polar  dipole  array  is  one  in  which  the  position  of  !he 

potendal  electrodes.  P,  and  P2,  in  Figure  A-2  remains  unchanged  but 

"L"  and  t^ed  t ' ^ “ PlaCCd  °"  ^ Vertical  ^averse  arm  of  the 
L and  the  distance  between  C,  and  C2  is  allowed  to  approach  zero 

Compavisons  of  these  ideal  array  configurations  have  been  analysed  by 


G.  V.  Kel.^.j  Dipole  Method  for  Deep  Resistivity  Studies  " 
Geophysics,  3^,  1088-1104,  1966. 
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probing  depth  sent!  tlv y has  * ^per 

*■  o'  an  epuat^wLpoie  a„r,  Ha”ay  "°'  25  «'“* 

resistivity  response  than  the  equatorial  dtnoleY  * C!"8"  aPParent 

as  that  of  a polar  dipole  array.  ) P rray  but  n0t  as  strong 

which  the  ^al  etcterod;SLTed  “ d«»“*d  ^ is  one  in 

shown  in  Figure  La  "T”  ?"en.ted  at  Points  ^ * and  P2 ' as 

arm  of  the  "L".  In  this  case^the'ellctr o^^  °rientati°n  alon8  the  upper 
current  distribution  produced  bv  the  so  ^ 15  resPonsive  to  the 

assumed  homogeneous  earth  ™ ^ electrod«.  C],  under 

Pi'  and  P,'  on  a ' I!"  co"^t1ons  because  of  the  location  of  both 
surround  C]  if  Q were  atTnT  bemis Phe rical  shell  which  would 

otherwise  rL  pons”!  to  the  ctr^t  iTTr  ^ ^ 

sink  electrode  and  can  be  used  asp  ri  U 10n  Produced  by  the  current 
of  the  current  sink  electrode  on  th  meaSU*‘e  of  the  perturbing  influence 
source  electrode  as  observed"  the  **  ^ 

procedure  is  useful  in  transforming  the  L shated  ThiS 

readings  to  those  obtainable  with  th*  s'  u , P d array  aPParent  resistivity 
this  method  the  infinite  distance  s'nk  T t l0Cated  ^ 

transformed  resistivity  will  corrps  *6Ctr°de  can  be  avoided  and  the 

in  Schlumberger  albumstOof  theoreti  "l  ° * 6 resistivity  data  available 
subsurface  cases.  theoretical  resistivity  curves  for  various 

L - shaped y°u s^ng d|  * C U S S 1 on  * the  more  generalized  use  of  the 
as  described  CE  “‘T*  * ^ Pa- 

components  of  app  ^ “d  "y" 
the  resultant  "total"  apparent  r*«icH  V . °f  these  components  and 
anisotropy  of  the  subsurface  earth  mlJer^s  Tnd  inf°rmati°n  °n  the 

array  can  bT  oplraUda^er^^a  ’ ‘”e  LH:ShaPed 

array- By  combinine  b°“-array 

advantages  In  certain  survey  app.i^I “a£”7' 

Geophysique,  ^ 

* £jXPt°r.  Geophysisists , 1963. 
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be  gained  by  this  method  are  in  eliminating  the  need  for  the  infinite  dis- 
tance current  sink  electrode  and  in  providing  interpretive  information 
on  the  anisotropic  resistiv.ly  signatures  of  elongated  solution  cavity 
structures  or  man-made  underground  mining  excavations  of  interest. 

V.  A SWITCHED-ELECTRODE  EQUATORIAL  DIPOLE  ARRAY 

In  a previous  investigation  pertaining  to  the  detection  of  shallow 
man-made  tunnels,  Southwest  Research  Institute  tested  an  unconvention- 
al electrical  resistivity  survey  method  intended  to  emphasize  the  re- 
sponse to  small  subsurface  tunnels  and  to  simplify  the  required  resis- 
tivity interpretation  requirements.  ^ The  work  on  this  method  was  a 
limited  effort  aimed  at:  ( 1 ) establishing  the  feasibility  of  the  method  for 
detecting  tunnels;  (2)  evaluating  its  performance  in  comparison  with  a 
standard  resistivity  method  (a  Wenner  array);  and  (3)  evaluating  the 
method  using  rolling- contact  metallic  electrodes. 

The  switched  - elect  rode  equatorial  dipole  array  was  first  con- 
ceived as  a vehicle  wheelbase  - size  electrode  array  whereby  the  neces- 
sary current  source,  potential  measuring  instrument,  switching  circuits , 
and  data  recording  and  output  display  equipment  would  be  carried  as  part 
of  a single  vehicle-mounted  survey  system.  The  electrodes  of  the  array 
would  be  mounted  as  part  of  the  vehicle  wheels  as  a means  for  continuous 
mobile  resistivity  surveyin'  over  suitable  off-road  terrain.  With  this 
method,  the  penetration  depth  of  the  current  distribution  patterns  would 
only  be  comparable  with  the  wheelbase  and/or  track  width  dimensions 
of  the  electrode  array;  however,  such  depths  were  considered  adequate 
for  the  intended  tunnel  detection  application.  The  exceptional  results 
obtained  with  the  simplified  test  version  of  this  concept  in  detecting 
typical  man-made  tunnels  in  the  San  Antoi>:o  vicinity  were  encouraging 
enough  to  suggest  the  use  of  this  method  on  a larger  scale  for  sub- 
surface cavity  detection  related  to  other  survey  applications. 

Figure  A-3  illustrates  the  general  electrode  configuration  re- 
quired in  the  switched-electrode  array.  On  one  half  of  the  switching 
cycle  (solid  lines),  electrodes  E)  and  E4  comprise  a current  dipole  of 
spacing,  a,  and  electrodes  E2  and  E3  comprise  a potential  measuring 
dipole  of  spacing,  a,  separated  at  a distance,  b,  from  the  current 
dipole.  On  the  alternate  half  of  the  switching  cycle  (dotted  lines), 
electrodes  Ej  ^ud  E2  comprise  a current  dipole  of  spacing,  b,  and 
electrodes  E3  and  E4  comprise  a potential  measuring  dipole  of  spacing, 
b,  separated  at  a distance,  a,  from  the  current  dipole.  The  spatial 
orientation  of  the  equatorial  dipole  array  pattern  established  on  the 
alternate  switching  half  cycle  is  rotated  90  degrees  from  that 
established  on  the  first  half  cycle. 


Fountain,  loc.  cit. 
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The  geometric 
operation  are: 


factors  for  the  two  half 


cycles  of  switched 


For  e2  = current  source  electrode 
E4  = current  sink  electrode 

2 _ potential  electrodes, 

E3 
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and  for 


= current  source  electrode 
= current  sink  electrode 

= potential  electrodes  j 


For  the  case  of  a square  electrode  pattern,  the  two  geometric  factors 
are  equal.  That  is,  for  a = b, 


Ked  = ked  = 


rr 


(A  - 2 0) 


H ..  *■  ThC  SWltchlnS  Process  was  incorporated  into  the  original  resis- 

nate  tL  Vf  7 ^°ncept  for  two  reasons.  The  first  was  to  elimi- 
of  ah  ! 7 J ca!culating  and  interpreting  the  survey  data  in  terms 
of  absolute  observed  apparent  resistivity  when  only  the  anomalous  vari- 
ations were  of  interest.  This  could  be  achieved  through  the  switching 
process  by  subtracting  the  resistivity  values  derived  on  each  half  cycle 
o operation  to  obtain  a detected  resistivity  difference.  In  this  manner, 
the  possible  wide  range  of  magnitudes  of  resistivity  encountered  in  dif- 
ferent soils  would  not  be  involved  in  the  survey  operator's  interpreta- 
tion process.  The  second  reason  for  switched -electrode  operation  was 
o provide  a change  of  array  orientation  as  the  survey  progressed  in 

asneV°  A ^ r.esistivity  responses  from  more  than  one  observational 
aspect.  As  found  in  the  field  tests,  the  wheelbase  electrode  array  ex- 
hibited significantly  different  apparent  resistivities  in  the  presence  of  a 
tunnel  target  for  each  half  of  the  switching  cycle  because  of  the  elon- 

inE  thenelUrf  H tUnnCl  rClated  t0  thC  array  pattern-  Thus,  byswitch- 
g the  electrode  orientation  progressively  along  the  traverse,  ^he 

c hances  of  missing  an  arbitrarily  oriented  tunnel  target  because  of  an 
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inappropriate  array  aspect  would  be  reduced.  In  addition  to  these  two 
advantages  of  the  switching  process,  it  was  also  found  in  the  course  of 
the  field  tests  that,  for  90-degree  electrode  array  switching  as  illus- 
trated in  Figure  A -3,  the  apparent  resistivity  anomaly  functions  obtain- 
ed on  each  half  cycle  were  complementary  in  form.  Thus,  when  the  dif- 
ference between  the  two  resistivities  was  obtained,  it  resulted  in  a 
stronger  net  anomaly  than  that  exhibited  by  either  half  cycle  alone. 

The  average  difference  resistivity  obtained  by  the  method  of 
electrode  switching  which  includes  subtracting  the  apparent  resistivities 
observed  on  each  half  cycle  and  filtering  the  cyclic  switching  frequency 
components  is 

A°a,b  = P'P  = kED  - k'ed  (A  -2 1 ) 


for  the  electrode  switching  arrangement  shown  in  Figure  A-3  and  for 
the  case  of  constant  current,  I,  flowing  in  each  switched  current  dipole 
pair.  Substituting  the  geometric  factors  from  Equations  (A -19)  and 
( A — 1 9 ) gives 


tt  </a2  + b^ 
I 


- ML,  - a A0 ' 

/a2  +b2  -b  y a^  + b2  -a 


(A -22) 


Because  of  the  difference  in  array  factors  the  net  resistivity  is 
not  zero  in  Equation  (A-23)  for  the  case  of  a homogeneous  earth  medi- 
um; however,  this  difference  is,  on  the  average,  a constant  value  which 
changes  relatively  slowly  along  the  traverse  and  could  be  nulled  out  by 
subtracting  a longer  term  average  of  the  difference  resistivity.  An  u 
if  desired.  Pa> b 

Another  advantage  offered  by  the  rectangular  electrode  array 
lies  in  the  fact  that  each  half  of  cycle  of  switched  operation  offers  a 
slightly  different  resistivity  penetration  depth.  Thus,  the  probing  depth 
of  the  composite  switched -electrode  array  is  optimized  for  a wider 
range  of  target  depths  than  either  one  of  the  switched  arrays  alone. 

For  the  case  of  a square  electrode  array  pattern,  the  average 
difference  resistivity  obtained  by  the  signal  processing  technique  des- 
cribed above  is,  for  a = b, 

Apa«a  = (y/2  ~ ) i (A0-A0).  (A-23) 
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Thus,  for  a square  switched-electrode  array  pattern,  the  observed  dif- 
ference resistivity  is  proportional  to  the  difference  between  the  poten- 
tials measured  on  each  half  of  the  switching  cycle.  This  difference 
clearly  approaches  zero  'or  the  condition  of  a homogeneous  earth  mate- 
rial containing  no  anomalies.  However,  when  an  arbitrarily  oriented 
anomaly  such  as  a cavity  is  approached  along  a survey  traverse  the 
difference  between  the  apparent  resistivity  responses  of  the  two’arrays 
will  reveal  the  presence  of  such  a target.  Moreover,  this  target  indica- 
tion will  be  further  enhanced  by  any  anisotropic  resistivity  effects  such 
as  those  exhibited  by  an  elongated  cavity  structure. 
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